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ABSTRACT 
This  dissertation  is  an  integrated  method  of  designing  sampled- 
data  feedback  control  systems  to  required  time  response.     This  method 
is  distinguished  from  the  classical  frequency  design  methods  in  two 
ways  -  first,   the  closed- loop  rather  than  the  open-loop  function  serves 
as  the  starting  point  of  the  method;  and,    second,   the  procedure  invol- 
ves nearly  direct  synthesis  rather  than  trial  and  error.     To  implement 
the  design  procedures,  a  set  of  correlation  theorems  are  derived  which 
permit  interpretation  of  the  closed-loop  system  time  responses  in 
terms  of  the  pole-zero  configuration  of  the  closed-loop  system  function. 
Further,   a  graphical  method  is  presented  which  permits  rapid  solu- 
tion of  the  open- loop  poles  when  the  closed-loop  poles  are  known.     The 
techniques  and  results  of  an  analog  computer  study  are  presented  which 
serve  as  an  experimental  verification  of  the  methods  developed  in  this 
dissertation. 
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CHAPTER    I 
INTRODUCTION 

1. 1  General  Background 

The  past  few  years,   since  the  end  of  World  War  II,  has  witnessed 
a  tremendous  increase  in  industrial  and  military  applications  of  feed- 
back control  systems.     As  a  result  of  this  activity,  the  field  of  feed- 
back control  systems  has  developed  and  specialized  in  many  direc- 
tions.    One  important  class  of  systems  that  has  recently  been  sub- 
jected to  vigorous  investigation  is  known  as  sampled-data  systems (  or 
sampling  servomechanisms).     In  this  class  of  systems  one  or  more  .of 
the  variables  consists  of  trains  of  pulses  or  sequences  of  values.     These 
pulses  are  generally  available  at  equally  spaced  intervals  of  time  and 
no  information  is  received  between  two  consecutive  pulses.     Conven- 
tional systems,  or  systems  in  which  the  variables  are  always  continu- 
ous functions  of  time,    may,   in  fact,   be  considered  as  a  special  limit- 
ing case  of  sampled-data  control  systems,   where  the  sampling  interval, 
or  period,   approaches  zero.     However,   the  present  body  of  published 
sampled-data  control  system  theory  is  most  useful  when  applied  to  sys- 
tems that  have  sampling  rates  low  compared  to  the  maximum  response 
frequency  of  the  system.     It  is  to  that  type  of  system  that  the  findings 
presented  in  this  manuscript  are  directly  applicable. 

One  of  the  earliest  applications  of  sampled-data  systems  was  a 
tracking  system  including  a  search  radar.     In  this  case  the  antenna, 
rotating  at  a  constant  angular  velocity,    results  in  information  about 
the  location  of  the  object  which  is  available  only  as  the  antenna  sweeps 
past  the  object.     Time  multiplexing  of  telemetering  data  is  another  com- 
mon source  of  sampled-data.     The  increasing  trend  toward  systems 
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containing  a  digital  computer  as  one  element  is  furnishing  further 
examples. 

In  the  aforementioned  applications  the  processes  require  that  the 
information  be  in  sample  form.  In  other  cases  sampling  is  purposely 
introduced  although  the  signal  information  is  available  as  a  continuous 
time  variable.  Generally  this  is  done  for  the  express  purpose  of 
increasing  sensitivity  or  for  more  economical  use  of  equipment.  The 
electromechanical  galvanometer  and  chopper  bar  is  a  classical  exam- 
ple of  such  a  device  (1). 

Among  the  early  contributors  to  the  study  of  these  systems  were 
Oldenbourg  and  Sartorius  (2).     They  applied  difference  equations  to  the 
analysis  of  sampled-data  control  systems.     The  solutions  of  these  equa- 
tions yielded  the  absolute  stability  of  the  system  as  well  as  the  response 
to  various  inputs.     This  approach  offered  the  same  difficulties  as  exper- 
ienced in  the  classical  methods  of  using  linear  differential  equations 
for  describing  continuous  signal  feedback  control  systems.     The  method 
is  cumbersome  and  does  not  permit  relating  the  performance  of  the 
systems  to  its  constants.     As  a  result  of  these  difficulties,  attention 
has  been  given  to  developing  transform  calculus  methods  similar  to 
those  available  for  continuous  systems. 

Through  the  independent  efforts  of  Hurewicz  (3),   Stibitz  and  Shan- 
non (4),  transform  methods  were  evolved  which  gave  input  and  output 

Ts 
relations  of  the  system  as  a  function  of  z,    where   z   is  defined  as   e 

(s   is  the  Laplace  complex  variable.)    This  actually  is  a  special  case 
of  the  Laplace  transform  in  which  the  transform  is  a  rational  function 
of  z    .     Because  of  this,   the  inverse  transform  yields  the  time  solu- 
tion only  at  the  sampling  instants. 
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Ragazzine  and  Zadeh  (1),   pursuing  the  developments  of  Hurewicz, 
formalized  a  transform  calculus  now  known  as  the  z-transform.     Trans- 
form pairs  similar  to  Laplace  transform  pairs  were  made  available. 
Transfer  functions  in  terms  of  the  variable    z   were  shown  to  be  deriva- 
ble in  straight  forward  manner.     Even  more  important,  they  demon- 
strated how  sampled-data  could  be  analyzed  by  z-transform  methods. 
This  article  laid  a  firm  foundation  for  future  frequency  analysis  inves- 
tigation.    Significant  among  the  many  recent  contributors  are  Jury  (5) 
(6)  (7)  (8)  and  Lago  (9)  (10).     They  showed  that  methods  and  concepts 
associated  with  the  analysis  of  continuous  systems  could  with  some 
modification  be  directly  applicable  to  sampled-data  control  systems 
using  the  z-transform.     Correlation  was  established  between  the  fre- 
quency response  and  the  transient  response  of  some  second  order  sys- 
tems.    This  information  was  shown  to  be  useable  for  approximate  design 
of  higher  order  systems  where  dominance  holds. 

Linvill  (11)  showed  that  the  transfer  function  in  a  sampled-data  sys- 
tem could  be  represented  in  the  frequency  domain  by  an  infinite  series. 
Fortunately,  only  two  or  three  terms  are  generally  needed  for  a  good 
approximation.     The  labor  involved  is  about  the  same  as  with  z-  trans- 
form techniques  but  it  has  the  advantage  of  being  easily  applied  by 
those  familiar  with  the  standard  methods  of  analyzing  continuous  vari- 
able systems. 

As  the  previous  discussion  indicated,   most  investigators  of  sam- 
pled-data control  systems  have  concerned  themselves  with  analysis. 
They  have  shown  that  either  the  z-transform  or  Laplace  transform 
methods  can  be  applied  to  the  open-loop  system  and  readily  yield  infor- 
mation as  to  the  response  of  the  closed-loop  system.     Although  the 
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labor  involved  is  considerable  greater  than  for  continuous  systems  of 
the  same  complexity,  both  methods  are  straight  forward.     That  the 
methods  are  applied  to  the  open-loop  system  is  significant,    since  in 
general  this  is  known  in  advance  by  the  designer.     Furthermore,  com- 
pensating networks,  employed  to  improve  the  closed-loop  responses, 
are  added  as  part  of  the  open-loop  system.     To  date,   design  procedures 
are  largely  based  on  analysis  and  trial  and  error  methods.     Predic- 
tion of  high  order  system  time  domain  performance  is  based  on  study 
of  second  order  systems.     In  these  respects  the  design  of  sampled- 
data  systems  is  similar  to  the  commonly  employed  methods  for  design- 
ing continuous  closed- loop  systems. 

In  continuous  system  design  considerable  use  is  made  of  the  abi- 
lity to  subdivide  the  system  into  transfer  function  blocks,   the  require- 
ment being  that  each  block  not  affect  the  operation  of  the  adjacent  blocks. 
This  requirement  is  satisfied  by  many  physical  devices  and  consequently 
the  open-loop  functions  of  continuous  servo  systems  can  readily  be    rep- 
resented by  tandem  transfer  functions.     An  analogous  situation  exists 
with  respect  to  sampled- data  control  systems.     However,   in  this  case 
the  transfer  function  must  be  applied  to  all  components  between  output 
and  input  of  samplers.     If  a  system  employs  only  one  sampler  then  one 
z-transform  function  must  be  used  to  describe  the  entire  open-loop 
system.     Whenever  a  new  network,   such  as  compensating  network,   not 
employing  a  sampler,   is  introduced  into  the  system  an  entirely  new 
transform  must  be  recalculated  for  the  complete  system.     Even  the 
slightest  change  in  a  parameter  value  of  an  existing  component  will 
require  the  determination  of  a  new  z-transform.     (Fortunately,   the 
introduction  of  a  sampler  into  an  otherwise  linear  system  does  not 
affect  linearity;  consequently,   gain  enters  into  the  z-transform  as  a 
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multiplying  factor.  )    Because  of  these  complications,  the  locus  shap- 
ing on  the  Nyquist  diagram  or  Bode  diagram  requires  prodigious  amounts 
of  labor  when  compared  with  the  amount  of  work  required  for  a  simi- 
lar analysis  of  continuous  systems.     If  this  complication  could  be  removed 
design  of  sampled-data  control  systems  would  be  no  more  difficult  than 
for  continuous  ones. 

This  investigator  spent  considerable  effort  studying  this  problem 
in  hopes  of  finding  some  design  methods  which  would  circumvent  or 
ameliorate  this  complication.     It  readily  became  apparent  that  the  solu- 
tion was  not  to  be  found  by  changes  or  additions  to  present  methods 
which  modify  the  open- loop  transfer  function.     Further  progress  could 
only  be  had  by  adopting  a  completely  new  approach. 

Since  working  with  the  open-loop  transfer  function  necessitates 
considerable  labor  and  complications,   it  was  thought  that  perhaps  the 
closed-loop  transfer  function  might  offer  opportunities  for  design  with- 
out the  proportional  increase  in  labor  experienced  when  this  approach 
is  used  in  continuous  systems  (12).     Since  system  performance  is  gen- 
erally  specified  in  terms  of  closed-loop  performance,  this  approach 
would  actually  have  the  advantage  of  being  a  design  method  based  on 
systhesis  as  contrasted  with  the  present  methods  based  on  analysis. 
Investigation  along  this  line  produced  a  direct  synthesis  method  which 
permits  the  designer  to  shape  the  time  response  of  the  closed- loop 
system.     Also,   the  procedure  indicates  and  helps  to  resolve  possible 
conflict  between  the  performance  specifications  and  the  plant  (fixed 
components  in  the  system).     When  applied  to  higher  order  systems, 
this  method  requires  less  time  than  those  methods  whose  starting  point 
is  the  plant  transfer  function.     The  detailed  results  of  this  investigation 
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constitutes  the  contents  of  this  dissertation. 

1.2   Discussion  of  the  Dissertation 

After  a  few  introductory  comments  the  subject  material  of  this  dis- 
sertation will  be  discussed  by  chapters.     The  figures  showing  system 
time  responses  were  determined  by  an  analog  computer.     These  com- 
puter curves  were  found  to  check  out  the  analytical  calculations  for 
values  at  the  sampling  instants  exceedingly  well,   the  error,   in  general, 
being  less  then  reading  errors  of  the  4  centimeter  paper.     So  far  as 
possible,  the  nomenclature  used  in  this  dissertation  is  that  recom- 
mended by  the  AIEE  Standard  Subcommittee  on  Terminology  and  Nomen- 
clature of  the  Feedback  Control  Committee  (13).     Because  of  the  nature 
of  this  dissertation  additional  terminology  was  required. 

Chapter  II  presents  various  aspects  of  current  sampled-data  theory 
which  will  be  required  for  subsequent  developments  in  this  disserta- 
tion.    Because  of  the  newness  of  the  field,  important  developments  are 
found  scattered  throughout  the  literature  and  it  would  not  be  expected 
that  the  average  reader  have  available  all  the  references  from  which 
the  material  of  this  chapter  is  drawn.     Hence,   Chapter  II  presents  this 
material  in  more  detail  than  would  be  normally  expected.     For  the  rea- 
der familiar  with  the  field,   this  chapter  could  be  scanned  or  omitted 
completely.     No  new  contribution  is  made  in  this  chapter  but  the  author's 
explanations  may  add  further  insight  into  the  theory. 

Chapter  III  is  the  major  contribution  of  this  dissertation.     A  philo- 
sophy of  sampled-data  system  design  is  developed  which  is  nearly  pure 
systhesis.     With  this  method  the  designer  is  able  to  take  time  domain 
specifications  and  proceed  in  a  logical  manner  to  the  final  system. 
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The  method  permits  including  specific  components  as  part  of  the  spe- 
cifications. 

The  procedure  is  quasi  graphical  leaning  strongly  on  root  locus 
techniques  (14).     This  will  find  favor  with  practicing  engineers  who 
generally  prefer  graphical  over  analytical  methods  since  they  are  less 
subject  to  mistakes  and  usually  provide  more  circumspection. 

The  procedure  may  be  described  by  the  following  five  steps: 

(1)  The  poles  and  zeros  of  the  closed-loop  function  are  deter- 
mined from  the  time  domain  specifications. 

(2)  The  poles  of  the  open- loop  function  are  found  by  graphical 
methods. 

(3)  The  open- loop  poles  and  zeros  are  shifted  to  include  the 
important  poles  and  zeros  of  the  specified  components. 

(4)  The  new  closed-loop  poles  and  zeros  are  determined  gra- 
phically and  checked  against  the  specifications. 

(5)  The  appropriate  compensation  function  is  determined. 
Steps    (3)   and    (4)   may  be  eliminated  if  the  designer  feels  that 

the  system  will  operate  satisfactorily  with  a  compensator  cancelling 
out  all  the  poles  and  zeros  of  the  specified  components. 

Chapter  IV  develops  the  necessary  techniques  to  implement  the 
design  philosophy  described  in  Chapter  III.     Step  (1)    requires  a  set  of 
correlation  theorems  which  permits  interpretation  of  the  closed-loop 
system  time  responses,    such  as  peak  output,   steady  state  error,   set- 
tling time,   etc.   in  terms  of  the  pole-zero  configuration  of  the  closed- 
loop  system  frequency  function.     Some  correlation  theorems  are  avail- 
able in  the  literature    (6)    (15).     However,  the  set  of  error  coefficient 
theorems  and  the  settling  instant  theorem  presented  in  this  chapter 
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were  derived  by  the  author. 

The  terms  closed-loop  root  locus  and  open-loop  root  locus  have 
been  introduced  to  distinguish  between  the  former  which  is  the  well 
known  "Evans  Technique"    (14)   and  the  latter  which  is  a  method  of 
finding  the  open-loop  poles  from  the  closed- loop  function.     In  this  dis- 
sertation the  open- loop  root  locus  was  a  necessary  development  to 
implement  step   (2).     This  technique,   originated  by  the  author,   repre- 
sents an  important  addition  to  the  present  root  locus  methods  for  sam- 
pled-data  control  systems.     Step   (3)   presents  no  difficulties  while  step 
(4)  utilizes  the  standard  root  locus  techniques.     Current  literature  has 
numerous  instances  where  the  closed- loop  root  locus  or  Evans  locus 
has  been  applied  to  sampled-data  control  systems   (6)   (12)   (15).     Only 
portions  of  these  references  which  are  directly  applicable  to  this  design 
method  are  included.     The  details  of  step   (5)  are  left  to  Chapter  V. 

Basic  pole- zero  configurations  are  derived  which  give  the  engin- 
eer a  starting  point  for  his  design.     It  is  shown  that  most  specifica- 
tions can  be  satisfied  by  a  basic  pole-zero  configuration  consisting  of 
one  pair  of  complex  dominant  poles  plus  several  additional  sets  of  real 
poles  and  zeros  which  must  be  carefully  located  according  to  the  rules 
developed  in  this  dissertation. 

Chapter  V   applies  the  technique  developed  in  Chapters    III  and 
IV  to  a  wide  variety  of  sampled-data  control  systems  which  included 
hold  networks  and  delays.     Although  not  an  exhaustive  list,  the  success 
of  the  design  methods  on  these  problems  indicates  its  general  applica- 
bility.    The  steps  involved  in  deriving  the  Laplace  transform  from  the 
z -transform  are  illustrated.     To  the  author's  knowledge  this  is  the 
first  time  that  this  process  has  been  presented  in  such  detail.     The 
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experimental  findings  verifying  the  methods  of  this  dissertation  are 
presented. 

The  general  theory  is  extended  to  types  of  sampled-data  control 
systems  which  do  not  enjoy  a  closed- loop  transfer  function  independent 
of  the  input  function.     This  is  a  further  complication  found  in  sampled- 
data  systems  not  having  a  direct  analog  in  continuous  linear  systems. 
These  systems  are  physically  characterized  by  the  omission  of  a  sam- 
pler in  the  error  channel.     Conventional  analysis  methods  readily  test 
the  absolute  stability.     However,   to  the  author's  knowledge  this  disser- 
tation is  the  first  to  include  design  procedure  for  these  very  special 
types  of  control  systems. 

Methods  for  handling  systems  with  multiple  inputs,   such  as  a 
disturbance  input,   are  also  discussed. 

Chapter  VI  presents  the  analog  computer  set-up  that  was  used 
to  verify  the  methods  in  this  thesis.  The  use  of  an  analog  computer, 
with  the  transfer  function  of  the  various  components  making  up  the  sys- 
tem directly  analoged,  permits  the  system  to  be  analyzed  without  the 
effects  of  nonlinearity  overshadowing  the  features  of  the  theory  to  be 
tested  in  the  experiment.  A  specially  designed  mechanical  sampler 
(16)   and  several  new  experimental  techniques  are  discussed. 


CHAPTER   II 
INTRODUCTION  TO  SAMPLED-DATA  THEORY 

2.  1      Operation  of  the  Sampler 

A  linear  sampled- data  feedback  system  can  be  represented  by  the 
introduction  of  a  switch  in  the  otherwise  continuous  linear  system,   as 
shown  in  Figure    2.  1.     This  switch,   referred  to  as  the  sampler,   is  the 
component  that  distinguishes  a  sampled-data  system  from  a  continu- 
ous system.     Figure   2.  1    is  misleading  regarding  the  operation  of  the 
sampler,   for  in  physical  systems  the  sampler  operates  in  accordance 
with  the  circuit  shown  in  Figure   2.  2   .     Here  the  switch,   SW,  is  alter- 
nately closed  on  position  a  for    T      seconds  and  closed  on  position  b 
for   T,     seconds.     The  travel  time  between  positions    a  and  b   is  con- 
sidered negligible  and  the  sampling  period   T   becomes 

T      =      T      +  Tb  (2.  1) 

In  this  dissertation  the  sampling  period  is  assumed  constant.     If  more 
than  one  sampler  is  present  in  the  system,   they  are  considered  as 
operating  with  the  same  period  and  synchronized. 

With  the  sampler  illustrated  in  Figure    2.2    ,  the  output  signal  is 
other  than  zero  only  for  the   T      seconds  each  period  the  switch  is  in 
position  a   .     Thus  the  operation  of  the  sampler  converts  a  continuous 
signal  e(t)   at  its  input  into  an  output  signal  of  pulses  with  width   T 
and  a  height  proportional  to  the  amplitude  of  e(t)   at  the  sampling  instant. 
Figure   2.  3   illustrates  the  input  and  output  signals  of  a  physical  sam- 
pler. 

In  the  systems  to  be  considered  in  this  dissertation,   the  time    T 
7  a 

is  much  smaller  than  both  T,     and  the  time  constants  of  the  systems. 
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Figure   2.1.       A  sampled- data  system 
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Figure   2.  2.       Physical  circuit  for  sampler 
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(a)    Continuous  signal  applied  to  the  input  of  the 
sampler 
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(b)    Pulse  train  at  output  of  the  sampler 
Figure   2.  3.       Input  and  output  signals  of  a  physical  sampler 
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Physically  this  latter  requirement  means  that  the  system  responds  a 
negligibly  small  amount  while  the  pulse  is  applied.     If  this  is  the  case, 
the  system  response  would  be  pie  same  to  an  impulse  with  strength 
equal  to  the  area  of  the  pulse.     Since  mathematically  the  impulse  is 
much  simpler  to  handle  than  the  pulse,   both  in  the  literature  and  in 
this  dissertation,  the  impulse  train  is  substituted  for  the  actual  value 
of  the  continuous  time  variable.     It  should  be  noted,   that  the  area  of 
the  impulse  is  considered  as  equal  to  the  value  of  the  continuous  vari- 
able at  the  instant  of  sampling  rather  than  corresponding  to  the  area 
of  the  pulse.     With  the  area  of  the  pulse  equal  to  the  value  of  the  con- 
tinuous variable  multiplied  by   T      ,  the  mathematical  representation 

a 

of  the  sampler  involves  a  gain  change.  Thus,  the  gain  K  of  the  open- 
loop  system  is  related  to  the  actual  gain  K  used  in  the  physical  sys- 
tem by: 

K 

r 

a 


Kq      =      -£-  (2.2) 


The  Figure   2.  4   may  yield  more  insight  to  the  nature  of  the  approx- 
imation involved  when  the  pulse  is  replaced  by  an  impulse  under  the  con- 
ditions stated  above.     It  is  clearly  shown  in  this  figure  that  the  fre- 
quency content  of  the  pulse  is  nearly  constant  over  the  extent  of  the 
bandwidth  of  the  open- loop  function.     Thus,   the  system  output  contains 
essentially  the  same  proportion  of  frequency  components  with  either 
the  pulse  or  impulse  input.     The  change  in  gain  is  illustrated  with  T 

ct 

less  than  one  second. 

It  is  helpful  for  analysis  to  consider  the  sampler  as  an  impulse 
generator  which  is  modulated  by  the  input  signal.     The  process  of  con- 
sidering sampling  as  modulation  of  an  impulse  train  is  shown  in  Fig- 
ure  2.  5    ,   wherein  the  values  of  the  impulses  are  represented  by  the 
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(a)      The  pulse  and  impulse  time  functions 


(b)    The  frequency  spectrums  of  the  impulse  (1)    ,   pulse 
(2)  and  bandwidth  of  the  open- loop  system   (3) 


Figure   2.4.     The  impulse  approximation  illustrated  in  terms 
of  the  Frequency  Domain. 
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(a)      Input  to  sampler 
i(t) 


/    - 


time 


r       ^r       sr 


time 


(b)      Impulse  train 


T  37        ST 


time 


(c)      Output  of  sampler 
Figure   2.  5   .     The  sampler  as  an  impulse  modulator 
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the  amplitude  of  the  lines  occurring  at  the  sampling  instants.  With 
input  and  output  denoted  as  e(t)  and  e*(t)  respectively,  the  output 
may  be  written  as 

e*(t)      =      e(t)i(t)  (2.3) 

where    i(t)    represents  a  train  of  unit  impulses 

i(t)      =     ]T     Uo  (t  "  nT)  (2'4) 

n=-  oo 

Equation  2.  3    may  be  written  equivalently  as 

e*(t)      =      £     e(nT)  uq  (t  -  nT)  (2.5) 

nao 

Here  the  negative  values  of    n    are  omitted  by  virtue  of  the  assump- 
tion that    e(t)    is  zero  for  all  negative  values  of  time.     The  Laplace 
transformation  of  both  sides  of  Equation   2.  5   yields 

oo 
E*(s)      *      £     e(nT)  e"  nTs  (2.6) 

n=o 

An  alternate  and  very  useful  expression  for      E*(s)      can  be 

obtained  by  the  use  of  contour  integration.     Equation  2.  3  represents 

the  sampled  time  function  as  the  product  of  two  time  functions  from 

which  the  Laplace  transform  can  be  obtained  by  a  process  of  complex 

convolution.     The  complex  convolution  integral  is  given  by  (17) 

c  +  j»« 

L  [fj  (t)  f2  (t)  ]    =    TJp      /      F^s  -  w)  F2  (w)  dw  (2.  7) 

c  -  j  oo 
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The  contour  of  the  integration  of  Equation   2.  7   is  a  path  from    c  -  j  oo 
to      c  +  j°°  along  an  analytic  strip  between  the  singularities  of 
F.   (8  -  w)      and  those  of      F^  (w)   . 

Letting  the  impulse  train  be  represented  by     f.   (t)      of  Equation 
2.7,       F.   (s  -  w)    becomes 


o-o 


F(s  -  w)   =   I(s  -  w)  =   Yj    e"  nt(8"w)  (2.8) 


n*o 

5 


This  infinite  summation  can  be  expressed  in  closed  form  as 

I(s-w)      =     1 (2.9) 

1   .   e-T(s-w) 

Then,  equating      e(t)      to     £y  (t)      yields 

F2  (w)      a     E  (w)  (2.10) 

With  the  indicated  substitutions  in  Equation   2.  7  the  Laplace  transform 
of     e*(t)     becomes 

re  +  j  oo 

L[e*(t)]     r       1-      /        E  (w)    1 dw  (2.11) 

ZWJ      J  ,         -  T(s-w) 

c  -  j  oo  l       e 

Evaluation  of  Equation   2,  11    can  be  effected  by  closing  the  poles 
I  (s-w)    in  the  right  half  of  the   w- plane  as  shown  inRgure  2.  6  .     When 
the  degree  of  the  denominator  of    E(s)    in    s    is  higher  than  that  of  the 
numerator  by  at  least  two,  the  integral  is  zero  along  the  infinite  semi- 
circle.    The  alternate  path  in  the  left  hand  plane  will  be  considered  in 
a  later  section  of  this  chapter. 

Applying  Cauchy's  integral  formula  gives 
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Figure   2.6.     Paths  of  integration  used  to  evaluate  the 
complex-convolution  integral. 
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Jrf 


E(w) 


1 


—    dw  =   -  /      Res 

.       -T(s-w)  f"-      , 

1-e  poles  of       |_ 

I(s-w) 


E(wn). 


.       -T(s-w   ) 
1-e  n  J 


(2.  12) 


The  poles  of  the    I(s  -  w)      are  the  zeros  of 


1    .    e-  T(s-w)      m      0 


(2.13) 


or 


l    ,  .    2  it  n 

U       =      s    +   j -— 

n  T 


(2.14) 


These  poles,   though  infinite  in  number,   are  all  simple  and  hence  the 


residues  are  given  by  Expression  2.  15 
n=  oo 


I 


Ejw) 


n-°°       d/dw[l-e-T(8"w)  ] 


(2.15) 


Since 


d/dw 


l.e-T(s-w) 


.  2irn 
<xy=s+j    ^ 


T   j  2  tr  n 
-  Te** 


(2.16) 


The  value  of  the  integral  represented  by  Equation  2.  12   is 


oo 


E*(s)    =    -L  ^     E(s  +  j 


2  ir  n 


ns-  oo 


(2.17) 


Defining  the  sampling  frequency  as   LJ     where 


u 


2    TT    /    T 


(2.18) 


Equation  2.  17   may  be  written  as 
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E*(s)      =     -L_     ^      E(s    +jn(Js)  (2.19)* 

T      n=-<x> 

It  is  seen  that  two  equivalent  expressions  are  available  for  analy- 
tical description  of  the  sampler,  as  demonstrated  by  Equation  2.  16 
and  Equation  2.  19.     The  first  equation,   which  expresses  the  Laplace 
transform  of  the  sampler  output  in  terms  of  the  values  of  the  input 
function  at  the  sampling  instants,  is  particularly  helpful  in  gaining 
insight  into  the  physical  operation  of  the  sampler  and  its  effect  on  sys- 
tem perf or mance.     The  second  equation  places  in  evidence  the  fre- 
quency effect  of  the  sampler  on  the  input  signal. 

Replacing    s    by    s  +  j  m(jj     ,   where    m    is  any  integer,   in  Equa- 

tion  2.  19  or  its  alternate  expression  Equation   2.  6,    results  in  an 

expression  which  is  identical  with     E*  (s)   (1)    (18)   .     That  is,      E*  (s) 

is  a  periodic  function  of    s    with  period      j  \jj     .     Expressed  symbo- 

lically 

E*  (s   +   j  m(J)    =    E*  (s)  (2.20) 

s 

This  indicates  that  the    s- plane  can  be  divided  into  strips  of  equal  width 
with  centers  at     n^       as  shown  in  Figure    2.  7  . 

The  spectrum  alteration  performed  by  the  sampler  is  shown  in 
Figure  2.  8.  From  this  it  is  seen  that  the  sampler  contains  an  exact 
replica  of  the  input  as  well  as  an  infinite  number  of  additional  images 


♦Equation   2.  19  is  in  error  if    e(t)  4   0    at    t  =  Of  .     In  this  case  the 
Equation  2,  19  should  be  modified  to  read 


n=    oo 


E*  (s)    o    _L_     £     (s   +   jn   U j   +    i/2  e(0+) 


T 

ns  -     oo 
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Figure    2.7.       Period  strips  of     E*  (s) 
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(b)      Amplitude  spectrum  of      E*  (s) 


Figure   2.  8.       Input  and  output  amplitude  spectrum  of 
a  sampler. 
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centered  about  a  multiple  of  the  sampling  frequency.     These  additional 
components  are  referred  to  as  complementary  signals  (11)  .     Spaces 
between  images  would  indicate  that  it  is  possible  to  regain  the  exact 
signal  by  filtering.     Unfortunately,  a  satisfactory  filter  for  the  process 
would  introduce  too  much  time  delay.     Any  signal  with  frequencies 
greater  than  the  sampling  frequency  would  overlap  into  the  adjacent 
strip  with  the  result  that  the  information  would  be  lost.     This  conclu- 
sion is  equivalent  to  Shannon's  sampling  theorem  which  requires  that 
the  sampling  frequency  be  at  least  twice  the  highest  input  frequency. 

2.2      Elements  of  z-Transform  Theory 

If  the  path  of  integration  for  Equation   2.  11    had  been  closed  to 
the  left,  the  finite  number  of  poles  of      E(w)      would  have  been  enclosed 
and  there  would  be  a  finite  number  of  residues  to  evaluate.     This  alter- 
nate path  is  also  indicated  in  Figure   2.  6.     In  this  case 


L  [e*  (t)  ]    =    V  Rei 


poles  of 
E(w) 


E(w  J 


n         .  -T(s-w   ) 

1  -  e       *         n    J 


(2.21) 


The  use  of  Equation  2.21    may  be  illustrated  by  determining  the 
z- transform  of  the  sampled  sequence  of  a  unit  ramp. 

e(t)  .=      t  (2.22) 

i  - 


The  addition  is  a  result  of  the  difference  in  definition  of  the  sam- 
pling process  in  Equation   2.  6   which  implies  a  full  sampling  at    t*0 
and  Equation  2.  11  which  assumes  the  impules  to  be  an  even  function 
sampling  only  one- half  of  the  value    e(o)  .     In  feedback  control  sys- 
tems Equation  2.19   is  normally  adequate. 
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Laplace  transform  is  then 

E(s) 

B 

1 

2 

(2.23) 

s" 

E(s)    is  seen  to  have  a  second  order  pole  at  origin.     Performing  the 
operations  indicated  by  Equation  2.21    yields 

E*  (s)      =     !-§ ^  (2.24) 

(e         -  1) 

An  alternate  method  of  deriving  Equation  2.  24   employs  Equation  2.  6 
as  a  starting  point.     Hence 

e*(t)    =    TuQ(t  -  T)  +  2TuQ(t  -  2T)  +  3TuQ(t  -  3T)  +  .  .  .  (2.  25) 

The  Laplace  transform  of  this  series  of  samples  is 


E(s)  =  Te"Ts  +  2Te'2Ts  +  3Te"3Ts  +  ...       (2.26) 


Multiplication  of  both  sides  of  the  above  equation  by  yields 


■•JS)      .     _!_   ♦   _L-      ♦     -L_      +...  ,2.27) 

_     Ts  2Ts  3Ts  4Ts 

T  e  e  e  e 

sT 
Multiplication  of  both  sides  of  the  above  equation  by        de  and 

integration  gives 

SliSL    deaT      .     -2 L-      -      -i—  ...    +C       (2.28) 

TesT  eTs  e2Ts  e3Ts 


where      C      is  the  constant  of  integration.     The  series  on  the  right  of 
the  above  equation  may  be  rewritten  in  closed  form. 
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JEli*)      desT     = ! +     C  (2.29) 

TesT  esT-l 

sT 
Differentiation  of  Equation  2.29  with  respect  to      e  yields 

T  e8T 
E*  (s)      =     — (2.30 

/   sT       1Y2 
(e        -   1) 

sT 
which  is  identical  to  Equation  2.24.     The  substitution     z    =    e  into 

the  expressions  gives  rise  to  the  z-transform.     The  above  expression 

in  z-transform  notation  becomes 

E(z)      «  T  Z  (2.31) 

(z  -  l)2 

i 

By  these  methods  it  is  possible  to  derive  a  table  of  z-transforms 
(5).     Table    2.1  is  a  short  list  of  z-transforms. 

The  inversion  of  a  z-transform  function  to  obtain  the  time  res- 
ponse of  the  function  at  sampling  instants  can  be  effected  by  any   one 
of  three  ways: 

(1)    Evaluation  of  the  inversion  integral    (19) 


*-£ 


-£r-<$)    C(z)zn"1dz  (2.32) 

where       '  is  a  path  of  integration  enclosing  all  the  singular 
points  of     C(z) 

(2)  Partial  fraction  expansion  of  the  function  so  that  each  term 
can  be  identified  with  the  corresponding  time  function  of 
Table    2.  1. 

(3)  Series  expansion  in  descending  powers  of    z    in  which  the 
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Table   2.  1.     Short  List  of   z-  Transforms 


Time  funtt)t>n 

Lap  face 
Transform 

2-Traris'torn) 

/ 

U4t) 

1 

1 

^ 

<A.(iJ 

i 

2 
2-1 

3 

t 

t 

Tz 

H* 

erQt 

i 

S  +  CL 

B 

Z-£'qT 

*    a  may  be  complex 
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coefficient  of      1/z        corresponds  to  the  value  of  the  time 
function  at  the  nth  sampling  instant. 

2.3      z-Transfcr  Function 

Of  primary  importance  in  servo  analysis  are  the  input- output 
relations  of  networks.     This  relation  for  sampled  data  systems  will 
now  be  established  using  the  simple  systems  illustrated  in  Figure   2.9. 
which  shows  a  linear  network      W      pulsed  by  a  sampler.     In  addition 
to  the  desired  continuous  output  time  signal  a  fictitious  sampler  is  pro- 
vided to  generate  a  sampled  output  signal.     Although  not  actually  part 
of  the  signal  flow  in  the  system,  consideration  of  the  sampled  output, 
rather  than  the  continuous  output,    simplified  the  analysis  as  will  be 
s  hown. 

It  is  evident  from  Figure    2.  9   that  the  Laplace  transform  is 
given  by 

C(z)      =      G(s)   E*  (s)  (2.33 

where      G(s)      is  the  Laplace  transform  of  network      W     and     E*  (s) 
is  the  sampled  input.     Applying  Equation  2.  19    yields  the  transform 
of  the  pulse  output  of  the  output  sampler. 

C*  (s)    =    -^r-  £     G(s  +  j  n   U3)  E*  (s  +  j  n   Ug)       (2.  34) 

Recalling  the  periodicity  of  pulsed-transfer  function  as  expressed  by 
Equation  2.20    ,   the  above  expression  may  be  rewritten 


C*  (s)      =     -IfY,     E*  (s)  G(s  +  jn  (J    ) 


(2.35) 
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Figure   2.9.       A  pulsed  network 
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Figure   2.  10.       A  system  with  a  sampler  between 
networks. 
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Since      E*  (s)      is  a  factor  in  each  term  of  the  above  series,  it  may  be 
placed  before  the  summation,   hence 


oo 


C*  (s) 


E*  (s) 


If-Ys   G(s  +  jn  U 


n=  -<X5 


(2.36) 


Denoting  the  bracketed  term  by     G*  (s) 

e>o 


G*  (s)      =      "T"  X!      G(s  +  Jn   UJs) 


(2.37) 


n=-oo 


Equation  2.  36  may  be  rewritten 


C*  (s)      =      E*  (s)  G*  (s) 


In  the  z- transform  notation 


C(z)      =      E(z)G(z) 


(2.38) 


(2.39) 


It  is  seen  that  the  z- transform  of  the  output  of  a  network  is  simply  the 
z-transform  of  the  input  multiplied  by  the  z-transfer  function.     Fur- 
thermore,  the  z-transfer  function  is  derived  in  the  same  manner  as 
the  transform  of  a  time  function. 

An  alternate  expression  for     G*(s)      in  terms  of    g(t)    ,  the 
impulse  response  of    W     ,  can  be  derived  by  noting  the  similarity  in 
form  between  Equation  2.  19   and  Equation  2.37   given  below  for  com- 
parison 

oo 

E*  (s)      «      -^r-      £      E(s  +  jn   U8) 
n«-  oo 

oo 


G* 


(s)      =     4-  G(s  +  jn  Ua) 


ns-  oo 
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Taking  the  Laplace  transform  of  Equation   2.  3    yields 

E*  (s)      *      L  [e(t)     K    i(t)  ]  (2.40) 

Since      G(s)      is  the  Laplace  transform  of    g(t)    it  follows  at  once  that 

G*  (s)      =      L[  g(t)   *  i(t)  ]  (2.41) 

The  over-all  transfer  function  of  two  networks  separated  by  sam- 
plers is  readily  seen  to  be  equal  to  the  product  of  the  z-transform  for 
each  network.     This  is  illustrated  in  Figure   2.  10   .     If  no  sampling 
device  separates  the  networks,   the  input  to  the  second  network  is  the 
continuous  function  of  time  present  at  the  output  of  the  first  network. 
Then  the  corresponding  pulsed  (starred)  transfer  function  will  be 


na 


c«o 


G*  (s)    =    "4^Z     G^S  +  j  n    ^Ja)  G2(s  +  j  n  Ue)  (2'42) 

na-oo 

Hence  the  over- all  transfer  function  must  be  calculated  by  considering 
the  system  as  a  whole.     The  accepted  notation   (1)  for  this  case  is 

G*  (s)      =      Gj  G*2  (s)  (2.43) 

or  in  terms  of  z-transform  notation 

G  (z)      =     GL  G2  (z)  (2.44) 

The  transfer  function  of  a  system  with  feedback,  such  as  shown 
in  Figure  2.  11  may  be  derived  in  the  following  manner.  The  expres- 
sion for  the  error  is 

E(s)    =    R(s)  -H(s)C(s)  (2.45) 

The  output  is  directly  related  to  the  sampled  error  by 
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Figure   2.11.     A  simple  feedback  sampled- data  system 
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C(s)    =    E*  (s)  G(s)  (2.46) 

Equation  2.46   substituted  into  Equation   2.45   yields 

E(s)    =    R(s)  -H(s)G(s)E*  (s)  (2.47) 

Taking  the  starred  transfer  function  of  both  sides 

E*  (s)    =    R*  (s)  -HG*  (s)E*  (s)  (2.48) 

Solving  for      E*  (s) 

E*  (s)    a    R*  (s)  (2.49) 

1    +  HG*  (s) 

The  output  is  then 

C*(s)    =    R*(s)G(s)  (2.50) 

1  +  HG*  (s) 

Taking  the  starred  transfer  function  of  both  sides 

C*(s)    =         R*(s)G*(s)  (2.51) 

1   +   HG*  (s) 

and  in  z-transform  notation 

C(Z)    =         R(z)  G(z)  (2.52) 

1    +  HG(z) 

Table   2.2   gives  the  transforms  of  the  output  for  a  number  of 
basic  servo  configurations.     Table  entries    2,    6  and   7    show  the  input 
may  not  always  appear  as  a  simple  multiplicative  factor  in  the  expres- 
sion for    C(z).     In  this  case  it  is  not  always  possible  to  derive  an 
expression  for  the  feedback  system  transfer  function  independent  of 
the  input.     As  discussed  in  the  introduction  this  creates  a  type  of  sys- 
tem whose  analog  is  not  found  in  linear  systems  operating  wholely  on 
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Table  2.  2.     Output  Transforms  for  Basic  Sampled-Data  Systems 
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continuous  data. 

2.4      Final  and  Initial  Value  Theorems,   Error  Coefficients  and  Mapping 

The  output  of  a  system      C(z)      may  be  expanded  in  partial  fractions 

as 

D   z  D.z 

C(z)    =    —2—     +    -^jr    +     Ig-T    +...  (2.53) 

z  -   1  1  2 

z-e  z-e 

The  inverse  transform  into  the  time  domain  is  readily  determined  by 
entry  4   in  Table   2.  1 

-d  T  -d2T 

c(t)    =    Dq      +      Dje  +     D2e      L        +  ...  (2.54) 

Assuming  that  the  real  parts  of  the  d's  are  positive  all  terms  of  Equa- 
tion 2.  54  go  to  zero  as  time  goes  to  infinity  with  the  exception  of      D     . 
Therefore,  the  final  value  of    c(t)    may  be  found  directly  from      C(z) 
as  follows 

lim(z  -  1)  C(z)    =    Urn  c(t)  (2.55) 

z   -^-1  t 


Expanding    C(z)    in  descending  powers  of    z    gives 


Hl  H2 


C(z)    =    H        +     -    +     —V1-    +    ...  (2.56) 

O  fc 

z  z 

The  inverse  transform  into  the  time  domain  is 

c(t)    =    H     +H.U(t-T)      +      H?U(t  -  2T)    +...  (2.57) 

O  1       O  £      O 

Hence,      H        ,   the  initial  value  of    c(t)      ,   maybe  found  directly  from 
C(z)      by 

34 


lim    C(z)      ■      lim    c(t) 
z  — ^oo  t  -*"0 


(2.58) 


The  error  expression  for  single  loop  systems  with  unity  feed- 
back is 

R(z) 


E(z)      = 


1  +  G(z) 


For  a  step  input 


Thus 


R(z)      = 


E(z)      = 


z  -   1 


(z-1)  [1  +G(z)  ] 


(2.59) 


(2.60) 


(2.61) 


Applying  the  final  value  theorem  the  steady  state  error      e  is 

s  s 


..  z-1 

e         =    lim    

ss  z 

t— ?oo 


lim 
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(2.  62) 


(z-1)  (1  +  G(z)  )  z^l     1  +  G(z) 


Then  the  error  coefficient  for  a  sampled-data  system  may  be  defined 
as 
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In  like  manner  it  can  be  shown 


K* 


and 


K* 


T" 


lim      (z-1)  G(z) 
z->l 


lim      (z-irG(z) 
z  —>l 


(2.63) 


(2.64) 


(2.65) 


(2.  66) 


It  was  explained  in  Section  2.  1    that  the  s-plane  could  be  consid- 
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ered  as  divided  into  strips.  The  conformal  mapping  of  the  first  strip 
in  the  s -plane  is  the  unit  circle  in  the  z- plane.  This  and  mapping  of 
other  important  lines  are  shown  in  Figure   2.12  . 

2.  5      Response  between  the  Sampling  Instants 

The  output  of  a  sampled- data  system  is  continuous  even  though 
its  input  is  sampled,  but  the  z-transform  only  contains  information  of 
the  time  response  at  the  sampling  instants.     This  section  presents  an 
extension  of  the  z-transform  called  the  modified  z-transform  which 
permits  determining  the  output  of  a  network  between  sampling  instants. 

Consider  the  system  in  Figure   2.  13    which  introduces  a  ficti- 
tious delay  between  the  output  of  the  network  and  the  sampler  at  the 
output.       A    represents  the  fraction  part  of  the  delay  to  the  total  sam- 
pling period  and  hence  has  a  positive  value  between  zero  and  one.     The 
effect  of  this  delay  on  the  signal  presented  to  the  sampler  is  illustrated 
in  Figure   2.  14  .     For  reasons  which  will  become  clear  later  it  is  aus- 
picious to  define  the  change  of  variable 

A      =      1    -    m  (2.67) 

Equation  2.67  implies  that      mT      is  a  pure  advance.     A  system  incor- 
porating this  pure  advance  before  the  sampler  and  a  delay  of  one  sam- 
pling period  following  the  sampler,  as  shown  in    Figure   2.  15,   is  seen 
to  have  the  same  output  as  the  system  in  Figure   2.  13.     The  impulse 
response  of  the  network  and  pure  advance  is      g(t   +  mT)      where    g(t) 
is  the  impulse  response  of  the  network  alone.     Employing  Equation 
2.  41    it  follows: 

Cn  (z)      =      L[g(t  +  mT)xi(t)]  (2.68) 
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Figure   2.  12.      Conformai  mapping  according  to 
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Figure   2.  13.       Samp  led -data  system  with  fictitious  delay,     AT 


Figure    2.  14.     The  output  compared  with  the  delay  output  of 
system  in  Figure    2.  1 1  . 
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Figure    2.  15.     Sampled-data  system  with  output  equivalent  to 
that  of  system  in  Figure   2.  13 
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Equation  2.68  may  be  evaluated  using  the  complex  convolution  integral 

in  the  following  form: 

c+joo 

Ca^     "     riTT  I      G(w)emwT    1 dw  (2.69) 

3J  ._  ,         -T(s-w) 

c-j  co  1  -  e  ' 

Evaluating  the  integral  in  the  left  hand  plane  provides  a  closed 
form  for      C*     (z)   .     It  should  be  noted  that  the  substitution  represented 
by  Equation  2.  67  assured  the  convergence  of  the  integral  on  the  infin- 
ite semi-circle. 

Because  output     C.(t)      is  identical  with     C*     (t)     but  shifted  in 
time  by  one  period,   it  follows  from  the  real  translation  theorem   (17) 

C*j   (s)     =      e"sT  C*q  (s)  (2.70) 

Varying  the  value  of    m    from  zero  to  unity,  the  entire  response 
between  the  sampling  instants  can  be  found.     The  z-transform  evolved 
from  this  technique  is  called  the  modified  z-transform.     With  this  nota- 
tion the  transfer  function  of  the  system  in  Figure    2.  13    ,   including  the 
linear  network  and  delay    AT      ,   is 


at    z'1    )  Rei 


G(z,m)      =    z         /         Residue 

pole  8  of 
G(s) 


~ .   .     msT 
G(s)  e 

sT      -1 
1  -  e         z 


(2.71) 


Table   2.  3   is  a  short  list  of  modified  z-transforms  for  various 
forms  of     G(s)  . 

2.  6     Analysis  and  Design  of  Sampled- Data  Control  System  by 

Conventional  Techniques 

The  absolute  stability  of  sampled-data  control  systems  can  be 
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Table    2.3.      Short  List  of  Modified   z- Transforms 


Laplace  Transform  Modified  z- Transform 
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determined  in  much  the  same  manner  as  is  done  with  continuous  sys- 
tems.    For  continuous  systems  a  very  powerful  technique  has  been 
devised  involving  the  construction  of  asymptotic  diagrams  in  which 
log  modulus  and  phase  are  plotted  against  log  frequency.     This  is 
known  as  the  Bode  diagram  and  is  primarily  useful  when  the  transfer 
function  is  expressed  as  a  product  of  simple  factors  and  the  variable 
increases  indefinitely  in  magnitude.     That  latter  requirement  is  not 
characteristic  of  the  z-transform,   where  for  this  diagram,     z    is 
varied  along  the  unit  circle. 

A  Nyquist  diagram  may  be  constructed  by  plotting  the  gain  and 
phase  of  the  open-loop  transfer  function  as    z    is  varied  along  the  unit 
circle.     Closure  of  the  coordinate      - 1  +  jO      indicated  instability.     Con- 
struction of  appropriate      M        circles    (20)   gives  an  indication  of  the 
relative  stability. 

As  an  example  of  analysis  and  design  using  the  Nyquist  diagram, 
consider  the  system  shown  in  Figure    2.  16.     The  sampling  period  is 
0.  5    seconds.     The  first  step  in  the  analysis  involves  determining  the 
open- loop  z-transfer  function,      G(z).     This  requires  a  partial  fraction 
expansion  of     G(s)      ,   converting  each  term  to  z-transform  with  the 
aid  of  Table  2.  1    and  finally  a  collection  of  these  terms  over  a  common 
denominator.     These  steps  are  illustrated  below. 

G(s)      *     L2 (2.72) 


s(s  +  3. 

2) 

3.125 

s 

- 

3.125 
s  +  3.  2 

3.125 
z  -  1 

-    - 

3.  125 
z  -  .2 

q(s)    =    ^i=r    -      -1^3^  (2.73) 

aw    =    44^-  -      4^-  U.74) 
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Digure    2.16.       Samp  led -data  system  to  be  compensated. 
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Analysis  on  the  Nyquist  plot,   Figure    2.  17   shows  the  system  to 
be  unstable.     Let  the  design  require  a  peak  overshoot  of   50%  in  the 
time  response  while  maintaining  the  same  low  frequency  gain.     As  a 
first  try,   it  will  be  assumed  that  this  requirement  could  be  met  with 
phase  lag  network  compensation 

Gcz<8>   "   -75  ;  t  i.5  <2-76> 

The  new  z -transform  open-loop  transfer  function  must  be  deter- 
mined as  shown  in  Equation  2.  77  through  2.80  . 

r>     i   \  7.  5(s  +  2)       .,        . 

Gl  (s)     =       s(s+  1.5)(s  +  3.2)  <2'77> 

r     <a\     -       3*125  ktI2  1.655 

Gl  (s)     *     i s  +  1.5       '    s  +  3.2  (2-78) 

n      .    .  3.  125z  1.47z  1.  655z  t0    __. 

Gl  (z)     3     TTT~    -      z-  .472    "   — T2-        <2- 79) 

r     t„\  1.  78z(z  -  .  41) .-    ftm 

Gl  (z)     "       (z-l)(z-.472)(z-.2)  <2'80) 

The  resulting  Nyquist  diagram,   Figure   2.17    ,  gives  the  locus 
touching  the      M      =    2.  5  circle.     This  indicates  a  time  response  over- 
shoot of  53%  .     Trying  a  second  network  compensator 

Gc2<s>  =  -5-f4r  <2-81> 

The     G(s)     and     G(z)      of  the  new  open-loop  system  are  then 

G2<s>    "      .[.'*%!  l\.l)  <2-82> 
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r     /    i      _      1. 57z(z  -  .  376)  n    „,. 

G2  (z)      "       (z  -  l)(z  -  .666)(z  -  .2)  (2*83) 

The  resulting  Nyquist  plot,    Figure    2.  17    ,   shows  the  locus  to  be 
outside  the      M        equal  two  circle,   indicating  a  peak  overshoot  of  less 
than  44%  .     Hence,   the  design  is  completed  as  far  as  this  procedure 
is  concerned.     The  first  step  in  calculation  of  the  time  response  at 
sampling  instants  is  to  form  the  output   z-transforms  as  given  in  entry 
4  of  Table  2.  2.     This  is  done  in  two  steps,   first  by  forming  the  closed 
loop  transfer  function  and  then  followed  by  multiplication  of  the  trans- 
form of  the  input.     The  final  result  is 

C(z)      =     1.  57  z    -     .59  z (2.84) 

z4   -    1.24z3    +   ,58z2    -    ,46z    +    .121 

The  first  few  values  of  the  time  response  at  the  sampling  instants 
are  placed  in  evidence  by  dividing  the  denominator  into  the  numerator 

as 

c(z)    »   -LH    +   -L|L  +   ^L  +   24L  + . . .  (2.  ss) 

z  z  z  z 

Although  the  actual  peak  overshoot  is  outside  of  specifications,  the 
Nyquist  procedure  has  given  a  system  which  at  least  has  a  reasonably 
good  time  response. 

The  previous  design  problem  has  illustrated  the  definite  trial  and 
error  nature  of  design  of  sampled-data  control  systems  by  conventional 
techniques.     For  higher  order  systems,   the  labor  may  be  prohibitive 
for  designing  to  fixed  specifications. 

Root  locus  methods  starting  with  open-loop  transfer  functions 
suffer  from  the  same  difficulties  in  determining  the  effect  of  a  com- 
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pensation  network. 

Examination  of  the  open- loop  transfer  functions  given  by  Equa- 
tions   (2.  75)    (2.  80)    (2.  83)   permits  drawing  four  important  conclusions 
about   z-transfer  functions  derived  from  continuous  components. 

(1)  The    s -plane  poles  have  direct  images  in  the  z- plane 

(2)  The    s-plane  zeros  do  not  have  direct  images  in  the    z-plane 

(3)  Additional   z-plane  zeros  are  introduced  when  taking  the 
z  -transform 

(4)  In  a  z-transform  the  numerator  is  of  one  degree  less  than 
the  denominator  if  the  function  has  a  zero  value  at  zero 
time  and  a  nonzero  value  at  the  first  sampling  instant. 
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Figure   3.  1.      A  compensated  sampled- data  control  system 


Figure   3.  2.       A  sampled-data  control  system  employing  a 
sampler  with  the  compensator 
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Figure    3.  1.       A  compensated  sampled- data  control  system 


Figure   3.  2.       A  sampled-data  control  system  employing  a 
sampler  with  the  compensator 
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culties  the  classical  approach  has  been  to  work  directly  with  the  open- 
loop  function.     A  compensator  is  chosen  to  reshape  the  Bode  diagram 
according  to  certain  tenuous  relationships  between  the  frequency  and 
time  domain.     For  continuous  control  systems  the  uncertainty  in  the 
final  performance  of  the  closed-loop  system  is  compensated  for  by  the 
ease  with  which  the  method  can  be  applied. 

The  design  example  in  Section  2.  6   showed  that  this  ease  of  pro- 
cedure for  design  based  on  the  open- loop  function  does  not  apply  to 
sampled  systems  unless  a  sampler  is  included  between  the  compensa- 
tor and  the  plant.     Such  a  case  is  shown  in  Figure   3.  2.     It  is  these 
difficulties  which  lead  to  the  philosophy  of  a  design  procedure  with  the 
closed- loop  functions  as  a  starting  point.     The  utilization  of  this  philo- 
sophy requires  two  new  developments: 

(1)  Method  of  determining  the  closed-loop  transfer  function 
from  the  specifications 

(2)  Method  of  determining  the  corresponding  open-loop  trans- 
fer function. 

This  dissertation  evolves  these  two  developments  making  possible  a 
logical  design  procedure  enabling  the  designer  to  take  the  specifica- 
tions and  proceed  in  a  direct  path  to  the  final  system. 

3.2      Synthesis  through  the  Closed- Loop  Poles-Zeros  Configuration 

The  details  of  the  synthesis  procedure  are  shown  by  the  flow 
chart  of  Figure    3.3.     The  various  blocks  represent  steps  in  the  pro- 
cedure.    The  letters  beside  each  arrow  identify  the  operation 
required  to  advance  to  the  next  block. 

The  chart  is  headed  with  the  time  domain  specifications.     By 
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use  of  "correlation  theorems"  it  is  possible  to  interpret  these  speci- 
fications in  terms  of  closed-loop  pole-zero  configuration  on  the  z-plane. 
From  the  pole-zero  configuration  the  "open  loop  root  locus"  yields  the 
corresponding  open- loop  pole-zero  configuration.     At  this  juncture  the 
designer  can  choose  either  of  two  concluding  procedures. 

(1)  Proceed  directly  to  determination  of  the  compensator 

(2)  Attempt  to  conciliate  the  differences  between  the  poles 
and  zeros  of  the  open- loop  function  and  those  of  the 

plant  so  that  at  least  the  major  poles  and  zeros  of  the 

> 

plant  are  included  in  the  final  open- loop  function. 

If  the  first  method  is  chosen  the  procedure  is  straight  forward. 

The  desired  open-loop  z-transform  is  converted  into  its  equivalent 

Laplace  transfer  function.     The  transfer  function  of  the  compensator 

.  may  now  be  found  by  simple  division 

G.  (s) 

Gr  (s)      =     — (3.1) 

Gp(s) 

where      G,  (s)      is  the  open- loop  or  "desired"  transfer  function  cor- 
responding to  the  selected  closed-loop  function.     The  actual  compen- 
sating network  may  be  determined  by  standard  network  synthesis 
procedure    (12) . 

A  study  of  Equation   3.  1    reveals  that  the  compensator  will,  in 
general,  cancel  out  each  pole  and  zero  of  the  plant  as  well  as  add  the 
poles  and  zeros  of  the  desired  function.     In  practice  it  is  impossible 
to  exactly  cancel  the  poles  and  zeros  of  physical  components  such  as 
electric  and  hydraulic  motors  since  their  precise  characteristics  are 
never  known  over  the  full  range  of  operation.     Fortunately,  complete 
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Figure    3.3.       Synthesis  Flow  Chart 
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cancellation  is  not  necessary.     System  performance  is  not  significantly 
changed  if  the  pole  and  its  cancelling  zero  only  lie  near  each  other. 
This  situation  gives  rise  to  a  pole  and  zero  placed  close  together  in  the 
closed-loop  function  and  the  resulting  residue  of  the  transfer  function 
in  the  pole  is  negligible. 

A  more  serious  objection  rests  on  the  argument  that  it  is  not  pos- 
sible to  cancel  out  equipment  time  constants  and  expect  the  equipment 
to  function  in  a  linear  mode  over  a  large  range  of  input  values.     It  is 
unreasonable  to  expect  that  cancelling  out  a  motor  time  constant  and 
replacing  it  with  a  much  smaller  one  will  result  in  any  increase  in  the 
actual  maximum  obtainable  acceleration  or  that  even  the  motor  and 
gear  combination  will  be  able  to  follow  a  ramp  which  exceeds  its  max- 
imum speed.     Hence  for  this  reason  it  is  strongly  urged  that  the 
designer  proceed  according  to  method   (2)  . 

In  method   (2)   the  designer  compares  the  poles  and  zeros  of  the 
desired  open- loop  with  that  of  the  plant.     By  shifting  the  location  of 
the  closed-loop  poles-zeros  it  is  usually  possible  to  include  the  major 
poles  and  zeros  of  the  plant  into  the  desired  open- loop  transfer  func- 
tion.    The  open  and  closed-loop  root  locus  furnish  a  graphical  method 
for  bringing  about  this  correspondence.     The  impossibility  of  incor- 
porating the  major  components  of  the  plant  into  the  desired  open-loop 
function  does  not  imply  failure  of  this  design  method.     Actually,   this 
is  a  signal  to  the  designer  that  the  fixed  plant  is  not  capable  of  pro- 
viding the  necessary  performance  demanded  by  the  specifications. 
This  simply  means  that  the  plant  must  be  improved  or,   if  this  is  not 
practical,   the  specifications  must  be  changed  to  concede  the  realities 
of  the  situation.     Thus  this  method  increases  the  designer's  circum- 
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spection. 

If  the  compensator  can  be  cascaded  to  the  plant  through  a  sam- 
pler as  shown  in  Figure  3.  2  ,  the  preceding  steps  are  still  taken  and 
the  compensator  can  be  found  directly  from  the  z- transform  of  the 

system. 

G     (z) 
G     (z)      =     — S_  (3.2) 

Gp  (z) 


3.3      Time  Response  and  Basic  Pole-Zero  Configuration 

The  critical  step  in  this  procedure  involves  choosing  a  pole- 
zero  configuration  to  satisfy  the  time  requirements.     The  typical 
time  requirements  for  a  transient  response  to  a  step  is  shown  in  Fig- 
ure  3.4.     The  steady  state  error  to  a  ramp  and  parabolic  input  are 
shown  in  Figure    3.  5.     Numerous  pole-zero  arrangements  could  be 
found  to  satisfy  a  particular  set  of  time  specifications.     However,  the 
correlations  theorems  to  be  formulated  in  the  next  chapter  are  based 
on  a  dominant  pair  of  complex  poles  which  together  with  a  real  zero 
will  fairly  well  control  the  transient  response.     A  pole-zero  spaced 
close  together  and  near  the  point      1,0      in  the  z-plane  can  provide 
control  over  the  steady  state-error.     With  this,  the  basic  pole-zero 
configuration  in  the  z-plane  is  that  shown  in  Figure   3.  6  . 

It  is  not  to  be  inferred  from  Figure    3.6   that  all  specifications 
can  be  discharged  by  third  order  systems.     However,   in  many  cases 
the  synthesis  can  be  handled  initially  in  this  way.     As  a  last  step  in 
the  design  additional  poles  or  zeros  may  have  to  be  added  if  the  com- 
pensator is  to  be  physically  realized.     (This  is  discussed  in  more 
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Figure    3.4.      Typical  transient  response  to  a  unit  input 
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Figure   3.  5.       Steady  state  errors  to  ramp  or  parabolic  inputs 
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Figure   3.  6,       Basic  closed-loop  pole-zero  configuration 
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Figure   3.7.       Typical  closed-loop  pole-zero  configuration 
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detail  in  Chapter   V.  )    If  the  poles  and  zeros  have  large  real  parts  in 
the  s-plane  they  will  give  rise  to  z-plane  poles  and  zeros  near  the  ori- 
gin and  have  little  effect  on  the  previously  determined  results.     Figure 
3.  7   shows  a  typical  pole- zero  configuration. 

It  should  be  kept  in  mind  that  this  control  of  time  response 
extends  only  to  the  sampling  instants.     This  generally  is  satisfactory. 
However,  a  complete  description  of  the  output  between  sampling  instants 
may  be  secured  by  employing  the  modified  z-transform  technique  to  the 
final  system. 

Plants  with  poles  and  zeros  in  the  right  hand  half  of  the  s-plane 
require  special  consideration.     The  poles  will  have  images  outside  the 
unit  circle  in  the    z-plane  and  so  will  the  zeros  in  most  cases.     Rather 
than  attempt  to  cancel  these  poles  and  zeros  they  should  be  included 
as  part  of  the  final  open- loop  transfer  function. 


/ 
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CHAPTER    IV 
SYNTHESIS  TECHNIQUES 


4.  1      Introduction 

This  chapter  develops  techniques  to  implement  the  design  phil- 
osophy described  in  Chapter    III.     First,    simple  relations  between  the 
pole-zero  configuration  and  the  time  response  are  developed.     The 
desired  response,   as  shown  in  Figure    3.4  and  Figure    3.  5    ,   maybe 
specified  by  the  following  five  quantities: 

(1)  The  peak  sampling  instant,     n     ,   the  sampling  instant  with 
the  maximum  overshoot  for  a  step  input. 

(2)  The  peak  sampling  overshoot,      M    ,  the  amount  of  the 
first  overshoot  at  the  sampling  instant    n      for  a  step  input. 

(3)  The  settling  sampling  instant,    n       ,  the  sampling  instant 

s 

that  the  value  of  the  output  reaches  within   5%  of  its  final 
value  for  a  step  input. 

(4)  The  number  of  oscillation,       N   ,   for  the  interval  up  to  the 
settling  instant  for  a  step  input. 

(5)  The  steady- state  error  constant      K*      ,      K*  and     K* 

7  P  v  a 

for  a  unit  step,   unit  ramp  and  unit  parabolic  input  respec- 
tively. 
The  first  four  quantities  describe  the  transient  response  while  the  last 
gives  the  steady- state  error. 

Basic  pole- zero  configurations  such  as  shown  in  Figure    3.  7  are 
next  derived  giving  the  engineer  a  starting  point  for  his  design.     Fin- 
ally rules  for  applying  the  open- loop  and  the  closed- loop  root  locus 

are  presented. 

58 


4.2      Correlation  Theorems 

For  the  sampled- data  feedback  control  system  shown  in  Figure 
4.  1    the  closed-loop  system  function     W(z)    is 

gWt     »      K°£l  (4.1) 


w  (z)     "      1  +  HGU)    *         "~FTO 


where 


w 


Q(z)      =  7T    (z  -  q,)  (4.2) 

k«l 


m 


P(z)      =       ill      <Z-Pi)  <4-3> 

The  values  of    q,       and    p,    may  be  either  real  or  complex,     It 
is  assumed  that  the  expressions  do  not  contain  multiple  roots  and  that 
the  order  of      P(z)      is  greater  than  the  order     Q(z)    .     These  two 
assumptions  are  true  for  most  physical  systems.     When  the  input  is 
a  unit  step,   the  output  at  the  sampling  instants  becomes 


c(nTV     ■•    Z"1 


W(z) 


rT 


(4.4) 


where      Z*        denotes  the  inverse  transform.     Thus  employing 


Equation  2.  32 


fz?'llP(z)     3nd*  (4.5) 


where      |       is  the  contour  of  the  integration  in  the  z-plane  that  encloses 
all  the  singularities  of  the  integrand  in  Equation  4.5.     Evaluation  of 
the  contour  integral  gives 
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Figure  4.  1.      A  sampled-data  control  system 
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Figure  4.2.      Important  characteristics  of  the  output 
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■I 


where 


z         =      u      +    j  v 
o  o         J     o 


P'(zo)      = 


d  P(z) 

d  z 


(4.11) 


z=z. 


It  is  important  to  realize  that  the  assumption  in  Equation   4.  10 
does  not  neglect  all  other  roots  but  neglects  all  other  modes  because 
these  modes  die  out  quickly.      The  accuracy  of  this  assumption  will  be 
discussed  later. 

In  general,  the  peak  of  the  output  will  not  occur  at  a  sampling 
instant.     If    n    is  permitted  to  be  noninteger,   then  a  value  of    n    could 
be  found,    so  that  at  one  sampling  period  of  time  later  the  output  will 
be  repeated.     This  condition  is  illustrated  by  Figure    2.2   and  can  be 
expressed  as 


c(n  +  1) 


c(n)      =      0 


(4.12) 


The  actual  value  of    n      is  the  upper  integral  value  of  the  actual 
n    obtained.     In  the  derivation  which  follows  it  will  be  assumed  that 
the  solution  of  Equation   4.  12    results  in  an  integer  value  for    n    . 

Substituting  Equation   4.12    into    4.8   and  solving  for    n       yields* 


n 


-    ang  Q(z    )   +   ang  P(z  J 


(4.13) 


*  The  intermediate  steps  leading  to  Equation   4.  1  3  are  quite  involved. 
The  more  important  steps  are  as  follows:    the  constant  term,   the 
coefficient  and  the    nth  power  of     |z     |       cancel  yielding 
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Then  the  peak  sampling  instant  can  be  expressed  as 


0     n 
o     p 


rr 


(sum  of  angles  from  the  zero  to 

the  predominate  pole)  -   (sum  of 

angles  from  the  other  poles  to 

the  predominate  poles)  (4.  14) 


cos 


n  0q    +   ang  Q(zq)    -    ang  (zq  -   1)    -    ang  P'(zq) 


z         cos 

o  ' 


(n+l)QQ  +  ang  Q(zq)  -  ang  (zq  -   1)  -  ang  P'(zq) 


Expanding  the  right  hand  side  of  the  above  equation  by  the  sum  of 
angles  trigonometric  identity. 


cos 


nM  ang  Q(zq)  -  ang  (zq  -   1)  -  ang  P'(zq) 


|zq|    |cos 


n  0     +  ang  Q(z    )  -  ang  (z     -   1)  -  ang  P(z    )  cos  0 


-    sin 


n  0q  +  ang  Q(zq)  -  ang  (zq  -  1)  -  ang  P'(zq) 


sin  0 


The  following  trigonometric  relationships  exist  for  the  values  above 


u 


cos  0 


cot 


sin  0 


ang  (zq  -   1) 


u     -    1 
o 


Substituting  these  into  the  preceding  equation  yields 
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Thus,   it  may  be  concluded  that  additional  poles  increase    n        and  that 

p 

a  zero  decreases      n 

P 

The  first  overshoot  is  obtained  by  substituting  Equation   4.  13 
into   4.  10   .     Hence 


M    = 


r/        M2,      2,1/2 
[(uo-l)  +vQ    ] 


K  Q(zQ) 


P'(zo) 


■J"* 


sin 


(ang  zq-1) 


(4.15) 


recognizing 


sin 


ang  (z-1) 


[(u   -l)2  +  v   2] 
L  x    o      '  o    J 


1/2 


(4.16) 


Equation   4.15    simplifies  to 


2  v 


M      = 


t(uo-D2+vo2] 


K  Q(zo) 


P'  (.„) 


I11 


(4.17) 


The  accuracy  of  the  assumption  that  all  modes  other  than  that 
due  to  the  constant  term  and  the  predominate  poles  may  be  explained 
readily.     From  the  complete  solution,   Equation   4.8,   consider  a     cer 


,n 


tain  mode      D(P.)     .       This  mode  will  decrease  to  a  small  percent  of 


Tan 


n  0Q  +  ang  Q(zq)  -  ang  P'(zq)  -  ang  (zq  -   1) 


Tan 


ang  (zq  -  1) 


From  which  Equation  4.  13   follows  immediately. 
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its  initial  value  when      n      =    3    ,     and  the  coefficient    D    is  small  if  it 

P 

is  not  too  close  to  the  other  roots.     In  the  event  that  a  certain  mode 

must  be  accounted  for  the  contribution  of  this  mode  at    n       can  be 

P 

evaluated  from  the  general  Equation  4.  8. 

Unless  the  actual  value  of    n    from  Equation  4.  13   is  equal  to  an 
integer  an  error  is  introduced  giving  an  optimistic  value  of    M.     This 
is  because  the  value  of  the  sine  term  in  Equation   4.  15,   for  the  non- 
integer    n    ,    is  less  than  the  sine  term  would  be  at    n      .     This  error 

P 

is  small  unless  the  sampling  period  is  a  large  proportion  of  the  period 

of  oscillation  and    n    is  near  midway  between  two  integer  values. 

For  a  zero  displacement  error  system  the  constant  term  must 

equal  unity  when  the  system  responds  to  a  unit  step  of  position.     Hence 

m 

rr 

ui  (i  -  P.) 

K      =      i_  (4.18) 

w 

7T  (i  -  %) 

Substituting  Equation   4.  1  8   into    4.  17   the  peak  sampled  overshoot  can 

be  expressed  as 

M      =         (Product  of  distances  from  all  poles  to  the  point 
of  1.  0),    excluding  distances  from  two  predomin- 
ate poles  to  the  point    1.  0) /(product  of  distances 
from  all  poles  to  predominate  pole   u     +  j  v       , 
excluding  the  distance  between  predominate 
poles)  x  (product  of  distances  from  zeros  to  the 
predominate  pole    u     +  j  v     )/(product  of  dis- 
tance from  all  zeros  to  the  point      1.  0)   |z    |   p 

(4.19) 

It  is  evident  from  Equation   4.  19   that  the  smaller    n      the 

P 

larger      M.       Thus  the  choice  of    n      and    M    requires  a  compromise. 

The  effect  of  real  poles  on      M      is  evident  from  Figure   4.3.       For 

b.   c 
position  one  near  the  origin     — »- is  nearly  unity  and  the  effect 
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z -plane 


Figure   4.  3.       The  effect  of  a  real  pole  on     M      and     n 


z -plane 


Figure   4.4.       The  effect  of  a  real  zero  on     M      and    n 
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t>2  C 
of  the  pole  is  slight.     However,   at  position  two,         -r- ^>    1 

will  increase.     Finally,    a  pole  in  position  three  near      c      , 

b3c 

<      1  M      is  reduced  while    n       is  increased.     Similar  con- 


b3  a  p 

elusions  may  be  deduced  when  the  poles  are  complex  . 

The  effect  of  a  real  zero  at  various  positions  is  evident  from 


Figure   4.4.     As  the  zero  approaches  the  point    1.0    the  ratio  becomes 

b  c 

increases  wniie  reducing    i 

P 
very  near  the  point      1.0      is  not  desirable  since      M      increases  tre 


D    C 

larger  and        -r- increases  while  reducing    n     .      A   zero  placed 


mendously.     For  most  cases,   it  has  been  found  that  one  zero  placed 
between  the  origin  and  a  position  somewhat  less  than  the  real  part  of 
the  predominate  pole  will  result  in  a  satisfactory  transient  response. 
This  corresponds  to  using  a  phase  lead  network  in  conventional  design 
methods.     Other  poles  and  zeros  near  the  origin  tend  to  cancel  each 
other  and  have  little  effect  on  the  time  response.     Thus,   a  pole- zero 
pattern,    such  as  shown  in  Figure    4.  5    consisting  of  a  single  real  zero 
and  a  pair  of  predominate  poles  should  serve  as  the  starting  point  in 
the  synthesis  of  transient  response. 

The  quantity      n        can  be  readily  determined  by  considering  the 
second  term  in  Equation   4.  10.       Since  the  cosine  term  is  always  less 
than  unity,   the  output  will  be  certainly  less  than    5%  of  its  final  value 
value  when  the  following  equation  is  satisfied. 

'       l1*6                           0.05  ..    ,n, 

z  =        a)      i (4.20) 


o  ■ 

2  K 


TR^T 


TTTT 


Taking  the  logarithm  of  both  sides  yields 
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z- plane 


Figure   4.  5.       A  basic  pole -zero  configuration  for  adjustment 
of  transient  response 
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1.3  -  Log 


K 


Q(zQ) 


n 


PT^0) 


Log    |z    | 


(4.21) 


A  response  with  no  output  at  zero  time  requires  the  coefficient 
of  the  second  term  in  Equation  4.  10  to  equal  the  reciprocal  of  the 
cosine  term  with  angle  of 


f  ang  Q(zQ)    -    ang  (zQ  -   1)    -    ang  P'(zq)  ]         (4.22) 


For  the  pole-zero  pattern  of  the  type  shown  in  Figure    4.  5,    the  angles 
are  those  shown  in  Figure   4.  6   and  generally  the  expression  4.  22   is 
about    180      and  the  coefficient  term  is  nearly  unity.     In  this  case  Equa- 
tion 4.  21    can  be  replaced  by  an  approximated  but  very  useful  result 


n 


1.3 


log 


(4.23) 


The  solution  of  Equation   4.23    generally  yields  a  noninteger  value  for 

n        and  the  actual  value  is  the  upper  integer  value  of  the      n        obtained, 
s  s 

The  locus  of      z        for  constant      n        are  circles  on  the    z-plane  with 
os  r 

centers  on  the  origin.      The  circles  are  shown  in  Figure   4.  7     and   can 
serve  as  starting  point  for  locating  the  predominate  poles. 

The  number  of  oscillations  may  be  calculated  approximately  as 
follows 

settling  time  (4.  24) 


or 


N    = 


N      = 


period  of  damped  oscillation 


n     0 
s     o 


2    TT 
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(4.25) 


z- plane 


Figure  4.  6.       Various  important  angles 
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z- plane 


Figure   4.7.       z       circles  of  constant    n 
6  o  * 
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For  control  systems  with  steady-state  error  to  a  unit  step  input, 
the  constant  term  in  Equation   4.  10  differs  from  unity  by  the  extent 
of  this  error.     The  value  of    K    may  still  be  calculated  since  the  amount 
of  the  error  is  given  in  the  specifications.     However,   in  most  practical 
cases,   the  error  is  negligibly  small  and  the  same  procedure  previously 
discussed  may  be  followed. 

The  error  constants  of  interest  are  the  position,   velocity  and 
acceleration  constants.     These  were  developed  in  section   2.4   in  terms 
of      G(z)      and  repeated  below  for  convenience. 

K*        =      lim      G(z) 
P  z-^1 

K*        =      -L-     lim      (z  ~    1)  G(z) 

Z— >1 

K*        =     —L,        lim      (z  -  l)2  G(z) 


a 


T  z-^1 


Introducing  a  change  of  variable 

6      =      z  -  1  (4.26) 

The  error  coefficients  become 


K*        =        lim      G(0)  (4.2  7) 

P  6  ->0 


K*        =     4r      lim     6   G(9)  (4.28) 

1  2 

K*        =     _L_      lim      e       G(0)  (4.29) 

a  T        6  -?0 


In  terms  of  the  variable      0      Equations   4.  2   and  4.  3   are 
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w 


Q(G)      =         Tf    (0    +    1   -  q,  )  (4.30) 

k=l  K 


m 


P(G)      =     7T  (o  +  i  -  p4)  <4-31> 

i=i  1 

From  Equation   4.  1     ,   the  open  loop  transfer  function  of  the  sys- 
tem shown  in  Figure    4.  1    with      H(s)      equal  to  unity  is 


w 
K   kTT      (6   +    1    -    qk) 
G<e)      =      _- (4.32) 

Tf   <e  +  i  -  P.)  -  k  rt    (e  +  i  -  qk) 

i=l  k=l 


Substitution  of  Equation   4.  32    into  Equation   4.  27   and  taking  the  limit 
yields 

w 

K       Tf       (1    -   q,) 

K*  =     ^Zl (4.33) 

P  m  w 

TT       (1  -  p.)   -      K     Tf    (1  -  q.) 

i=l  l  k=l  R 

Increasing  the  value  of    K    will  effect  an  increase  of      K*         .     A 
second  approach  would  introduce  a  pole  near  the  position     1.  0    and   a 
zero  between  it  and  the  origin. 

Incorporating  the  substitution 


a1       =       r-rj—  U.34) 


\    -    -r^kr  <4-35' 


into  Equations    4.30   and   4.31    yields 
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w 


Q(G)      =       Tf      (0      +     r-L  )  (4.36) 

k=l  bk 


m 


P(6)        =      Tt     (9      +      -±-    )  (4.37) 

1-1  ' 

Introducing  the  change  of  variable  and  the  substitution  above  in  Equa- 
tion 4.  18   gives 

w 

k-1         bk 

K  =         — £li £_  (4.38) 

m 

Thus  it  follows,   for  a  system  with  unit  feedback  and  no  steady- state 
error  to  a  position  input,   the  closed-loop  function  becomes 


w 
W(0)  =       _£zi 1 (4.39] 


TT        (a.  e  +  i) 
i=l 


The  open- loop  transfer  function  is  then 


w 


TT      (bk  e  +  i) 

G(0)        =  k=1 (4.40) 


m  w 


77     (a.  e  +  i)  -  Tf    (b  e  +  i) 

i=l  k=l 


Expanding  the  above  equation 
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w 


6' 


G(9) 


W     .2    vf     " 


+  •  •  •  ±6TTbk 


d 


m  w 


e 
+  J 


.rv 


m 


Multiplying  the  above  equation  by     W:  and  then  setting    ©    to  zero  ai 
indicated  by  Equation  4.28   yields  the  velocity  error  constant 


T  K* 


w 


(4.42) 


2>" 

=  1 


m 


w 


m 


If       /       b      -      )      a  is  zero  the  velocity  error  constant  is  infinite 


k=l  i=l 

and  the  acceleration  error  constant  becomes 


T2  K* 


w 


m 


L        k  ~  i 

c»  1  i=  1 


(4.43) 


In  terms  of  the  z-plane  poles  and  zeros 


1/T 


K* 


m 


i,f  Upi 


w 

v 

k=l 


^ 


(4.44) 
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2    /  T2 

K*         =        —L± (4.45) 

a  w  m 


k=l   \      k      i=l        pl 


The  denominator  of  the  above  equation  can  be  reduced  (large     K*    ) 
by  placing  a  zero  near  the  point   1.0   .     However,    this  will  have  the 
additional  effect  of  increasing    M    and    n  .     This  will  be  illustrated 

by  the  following  example.     Consider  a  simple  system  with  a  configura- 
tion such  as  shown  in  Figure    4.  5    . 


K  (z  -  qQ) 

W(z)     =         : (4.46) 

(z-u     +  j  v    )  (z  -  u     -jv) 
x  o       J     o'  x  o       J     o' 


Placing  these  poles  and  zeros  of  the  closed-loop  function  into 
Equation  4.46  ,  placing  the  expression  to  zero,  and  solving  for  q 
which  would  yield  an  infinite  velocity  error  constant  gives 

(1  -  u    )2   +  v 

qo        =         1 2 2 (4.47) 

2(1   -  uo) 

Assigning  values  of      u       =0.8,         v       =0.4      results  in  a      q        of 
°        °  o  o  ^o 

0.  5    .     This  location  for  the  zero  will  yield  an  infinite  velocity  constant 
but  the  overshoot,       M    ,    for  this  system  is     0.  6   .     If  the  zero  had 
been  placed  at  the  origin  the  value  of      M      would  be  a  more  reasonable 
0.39      but  the  value  of      T  K*        would  be  unity.     Thus  it  appears  that 
compensating  with  one  zero  does  not  give  independent  control  over 
peak  overshoot  and  velocity  error  constant. 
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An  arrangement  which  does  permit  adjustment  of      K*        without 
appreciable  effect  on  the  transient  response  makes  use  of  a  pole  and 
a  zero  placed  very  close  together  in  the  form  of  a  dipole  and  near  to 
the  point    1.0   as  shown  in  Figure   4.8    .      The  zero  at    q.    and  the  pole 
at    p        are  placed  sufficiently  close  together  that  the  residue  of  the 


r 
transfer  function  in  the  pole  at      p        is  made  negligible.     Although    p 

and      q.      are  nearly  equal,   there  can  be  a  significant  contribution  to 

the  denominator  in  Equation   4.  44  because  the  inverse  values 

-» and        -, are  large  numbers.. 

1  -  q  1  -  P 

Consideration  of  the  Equation   4.45  shows  that  a  dipole,  but  with 
the  zero  closer  to  the  origin,    may  be  used  to  raise  the  acceleration 
constant.     Furthermore,   additional  zeros  would  be  required  for  stabi- 
lization purposes  giving  rise  to  a  typical  pole- zero  configuration  shown 
in  Figure   4. 9    . 

4.  3      The  Required  Number  of  Poles  and  Ze  ros 

The  general  location  of  the  poles  and  zeros  to  achieve  various 
requirements  was  shown  in  the  previous  section.     However,  at  the  start 
of  the  synthesis  procedure  it  is  necessary  to  know  the  minimum  num- 
ber of  poles  and  zeros    that  are  needed  by  a  given  application.     Although 
mathematically  permissable,   a  practical  compensator  cannot  possess 
poles  at  infinity.     Consequently,   the  excess  of  poles  over  zeros  for  the 
compensated  or  desired  open-loop  function  must  be  at  least  equal  to 
the  excess  for  the  plant      G    (s)    .     Since  the  closed- loop  function  zeros 
are  the  same  as  those  of  the  open- loop  and  the  number  of  poles  are  the 
same  for  either  open  or  closed  loop  function,   it  can  be  said  that 
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z -plane 


Figure   4.  8. 


Pole- zero  configuration  for  large  velocity  error 
constant 


z- plane 


Figure   4.9. 


Pole- zero  configuration  for  large  acceleration 
error  constant 
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(number  of  poles  -  number  of  zeros) 

wu)  ^ 

(4.48) 

(number  of  poles  -  number  of  zeros)    ~    .    , 

Gp(z) 

It  would  be  of  considerable  advantage  to  avoid  the  lengthy  pro- 
cess of  taking  the  z-transform  from  the  Laplace  transform  in  order 
to  determine  the  number  of  zeros  of  the  plant.     In  general  this  can  be 
done  by  inspection  as  will  be  illustrated  by  considering  the  open- loop 
systems  shown  in  Figure    4.  10.      When  the  continuous  system  as 
shown  in  Figure    4.  10a  has  two  poles  in  excess  of  the  zeros,    output 
will  be  available  at  the  second  sampling,   i.  e.  ,     n  =  1.     Since  the  value 
of  the  output  can  be  determined  at  the  sampling  instant  by  dividing  the 
denominator  of      G(z)      into  the  numerator,   an  output  at  this  sampling 
instant  requires  that  the  number  of  z-plane  zero  be  only  one  less  than 
the  number  of  z-plane  poles.     The  number  of  zeros  equals  the  number 
of  poles  when  the  Laplace  transform  poles  exceed  the  zeros  by  only 
one.     Also,   in  these  cases,  a  zero  is  found  at  the  origin.     If  a  hold 
network   (1)   is  employed,   as  shown  in  Figure     4.  10b,   the  same  pro- 
cess of  reasoning  indicates  that  the  number  of  poles  exceeds  that  of 
zeros  by  one;  however,   a  zero  need  not  necessarily  be  at  the  origin. 
A  plant  including  a  delay,    such  as  Figure   4.  10c    ,   of  less  than  one 
sampling  period  will  give  rise  to  the  same  z^-plane  pole-zero  differ- 
ence as  discussed  for  Figure   4.  10a   except  the  zero  will  not  gener- 
ally be  found  at  the  origin.     For  every  complete  sampling  period  that 
the  delay  involves,    poles  are  placed  at  the  origin  without  any  additional 
zeros.     The  overall  transfer  function  for  Figure   4.  lOd    is  the  product 
G    i   (z)    GD2  (z)      and  the  difference  between  the  poles  and  the  zeros 
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(a)      The  plant  is  a  continuous  network 


cw  , 


*  ■» 


(b)      The  plant  includes  a  hold  network 


E&    /  eoa 


Ct2) 


a*)     . 


(c)      The  plant  contains  a  pure  delay 


Gp^S> 


Cte) 


CCS) 


(d)      Parts  of  the  plant  are  separated  by  samplers 


Figure  4.  10.       Various  types  of  plents 
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of  each  function  is  additive  in  the  final  transfer  function.     If  the  poles 
and  zeros  differ  by  one,     then  the  over-all  transfer  function  will  have 
an  excess  of  two  poles  over  the  zeros.     This  is  in  accordance  with  phy- 
sical expectation  since  the  second  network  does  not  receive  an  input 
until  the  second  sampling  instant. 

4.4      Open- Loop  Root  Locus 

The  method  of  this  dissertation  assumes  the  closed- loop  function 
to  be  known  leaving  the  open- loop  function  to  be  determined.     For  a 
system  such  as  Figure    4.  1      the  open-loop  function     G(z)      and  the 
closed- loop  function      W(z)      are  related  by 

Solving  for      G(z)      gives 

G(2)    ■     i  -Ww(i)  <4-50' 

Hence,   the  poles  or  roots  of      G(z)      are  the  solution  of  the  equation 

1     -    W(z)      =      0  (4.51) 

Both  functions  are  rational  polynomials  in    z     so  that 

w 
K    Tt      (z  -  qk) 

^--rg^-     -£ «"> 

TC   (z  -  P.) 

i=i  l 


ri    \  k  A(z) 

G(z)    =    -BTzT  = 


k  toTi  <z  -  ck) 


m 

TT     (z  -  P.) 

i=l 


(4.53) 
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Equation   4.  50   may  be  also  written  as 

kA(z)  K  Q(z) 

B(z)        =    P  (z)    -     k  Q  (z) 


(4.  54) 


The  gain  constant  are  the  same,   the  zeros  of  the  open- loop  function 
are  identical  with  the  zeros  of  the  closed-loop  function  and  the  poles 
of  the  open-loop  function  are  the  value  of    z    for  which 

W(z)      =      ej  2  h  *  (4.55) 

The  locus  of  roots  consist  of  all  points  on  the    z- plane  for  which  the 
plane  of  the  closed- loop  is      2  h  it      (h   may  be  any  integer  or  zero) 
that  is, 

ang    [   W(z)  ]    =    0  +  2  h  it  (4.  56) 

In  terms  of  Equation   4.  52   the  phase  of       vV(z)      is  given  by 

w  m 

0      =    2_     ang  (z  -  qk)      -        2_    ang  (z  -  p.)  (4.5  7) 

k=l  i=l 

The  precise  position  of  the  poles  on  the  root  locus  requires  that 

I   W(z)   |       =      1  (4.58) 


or  that 


w 


k  TT      | 


z   -  q 


k 


— *      1  (4.59) 

m 

Tf      h-pj 

i=l 


The  first  step  in  applying  this  graphical  root  finding  method  is 
to  determine  the  locus  of  roots.     An  approximate  root  locus  can 
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generally  be  found  very  quickly  and  simply  by  applying  some  simple 
rules  concerning  the  geometry  of  the  zero  phase  lines  in  the  z-plane. 
After  the  approximate  locus  has  been  obtained  the  exact  location  of  the 
roots  can  be  determined  by  use  of  Equation   4.  57   and   4.  59*  .     A  typi- 
cal open-loop  root  locus  is  shown  in  Figure    4.  11 

The  rules  and  their  proofs  for  the  approximate  construction  of 
the  open-loop  root  locus  are  now  presented. 

Rule    1.    "The  loci  are  continuous  curves  starting  at  a  pole  for 
zero  gain  and  terminating  at  a  zero  at  infinite  gain.  " 

Rule    2.     "The  loci  exist  on  a  part  of  the  real  axis  where  an  even 
number  of  poles  plus  zeros  or  none  are  found  to  the  right.  " 

This  may  be  easily  verified  by  consideration  of  Equation   4.  55 
and  Figure   4.11. 

Rule    3.      If  the  number  of  poles,   w,    exceeds  the  number  of  zeros, 

m    ,     there  are      w  -  m      branches  which  go  to  infinity.     Further,   the 

directions  of  the  infinite  asymptotes  with  respect  to  the  real  axes  are 

**: 

(4.60) 


w  -  m 


where      n   =   0    ,     1    ,    2    ,    3 

This  rule  is  proven  by  reference  to  Figure    4.  12a  and  b.     The 
scale  of  the  axes  in  Figure    4.  12b  is  increased  over  that  of  Figure  4.  12a 
reducing  the  location  of  the  poles  and  zeros  of      W(z)      to  almost  a    point, 
As  point      P      moves  further  away  from  the  origin  the  angles  subtended 
by  the  zeros  or  poles    are  essentially  the  same.     The  sign  of  the  angles 
from  the  zeros  is  apposite  to  the  sign  from  the  poles.     Since  the  angles 


*For  this  operation  a  device  called  the  Spirule  has  been  developed  by 
Spirule  Company,   9728  El  Venado,    vVhittier,   Calif. 
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z-plane 

x  closed- loop  pole 

0  closed- loop  zero 


Figure   4.  11.      A  typical  open- loop  root  locus 
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(a) 


(b) 
Figure   4.  12.       Proof  of  rule  three 
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from  the  excess  poles  are  all  the  same,   the  angle  of  each  pole  is 
2ir 


w  -  m 
prove  useful. 


if  the  sum  is  equal  to      2  it    .         The  following  table  will 


TABLE   4.  1 
Direction  of  Asymptotes 


w   • 

■  m 

D 

irect 

ions 

1 

0 

2 

-180° 

0 

3 

-120° 

0 

±120° 

4 

-90° 

0 

+  90° 

+  180° 

Rule    4.       "All  asymptotes  intersect  at  a  point      u        on  the  real 


axis.  " 


m  w 

Y  Y 

La  p.      -       La       q, 

i=l  1             k=l        K 


where  u        =       (4.  61) 

w  -  m 

The  proof  of  this   rule  can  be  demonstrated  by  consideration  of 
the  transfer  function  of      W(z)      written  in  the  form 


m  m- 1 

W(z)    =    K       z       *   alz  *   "'        am  (4.62) 

W  ,  W-l  u 

z       +   b.  z  +   .  .  .     +  b 


Dividing  the  numerator  into  the  denominator  to  give 

K 

W(z)    =        = (4.63) 

w-m        ,,  x     w-m-1    ,  \    .       / 

z  +   (bj   -  aj)  z  +  . .  . 
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The  locus  is  described  by  equating  the  denominator  to      -  K,      that  is 


w-m  /v.  \     w-m-1  v  ..    ,.. 

z  +    (bj  -  a1 )  z  +  ...=      -  K  (4.  64) 


By  the  theory  of  equations  the  sum  of  the  roots  of  Equation  4.  64   is 
b.    -    a.       .     When  the  magnitude  of    z    is  large  Equation  4.64  behaves 
as  a  polynomial  of      w  -  m      degree  with  the  sum  of  the  roots  given  by 
b.    -    a.       being  a  constant  which  is  independent  of  gain  as  long  as 
w  -  m      is  greater  than  two.     This  is  no  limitation  as  an  inspection  of 
Table   4.  1    will  show.     With     b.   -  a.       independent  of  gain,   the  sunn  of 
the  root  remains  at  one  point  for  any  value  of  gain.     Hence,   the  asymp- 
totic lines  radiate  from  the  point  on  the  real  axis  given  by 

b      -    a 

u        =     _J L  (4.65) 

e 

w    -    m 

Rule    5.        "The  point      u        at  which  the  locus  breaks  away  from 
the  real  axis  is  found  by  proper  equating  of  the  sum  of  reciprocal  of 
the  distances  between  poles  and  zeros  and  the  point    u      as  explained 
by  the  following:  " 

The  appreciation  of  this  rule  and  its  proof  is  demonstrated  by 
Figure   4.  13.       A  trial  point  is  assumed  at    u      and  moved  a  small 
vertical  distance      L  v      off  the  real  axis.     To  be  a  point  on  the  root 
locus  this  point  must  satisfy  the  phase  angle  criterian.     That  is 


-  e     +  e2   +  e3    +  e4    =    2  h  it  (4.66) 


All  the  angles  are  small  and  can  be  replaced  by  their  tangents,  hence 
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Ue  ^^ 


Figure   4.  13.       Proof  of  Rule  Five 


rc,     *-* 


Figure   4.  14.       Proof  of  Rule  Six 
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*  v  +    180°   -        A  v       +    180°    =    360°    (4.67) 


U      -   U.  U      -   U->  u,    -   u,  u.    -  u 

el  e         2  3         2  4         e 


which  becomes 


u     -  u.  u,  -  u  u.  -  u  u     -  u-, 

el  3e  4e  e2 


(4.68) 


The  solution  of  Equation  4.68   for  the  value  of    u       is  generally 
accomplished  by  trial  and  error.     For  the  case  where  complex  poles 
or  zeros  exist,   the  expression  needed  to  determine  the  break  away 
point  is  more  complicated.     However,    it  has  been  shown  (21)  that  unless 
the  complex  poles  or  zeros  are  close  to  the  point  of  emergence  their 
contribution  is  slight. 

Rule    6.     "The  angle  at  which  a  locus  emanates  from  a  complex 
pole  or  terminates  upon  a  complex  zero  is  found  by  summing  the  angles 
to  all  the  other  poles  and  zeros.      When  emanating  from  a  pole  the 
angles  from  the  poles  are  taken  negative.      When  terminating  on  a  zero 
the  angles  from  zeros  are  taken  negative." 

Consider  a  test  point      P        shown  in  Figure   4.  14   as  on  the  locus 
only  a  very  small  distance  removed  from  the  complex  pole.     Then  the 
angle      P        can  be  considered  the  angle  of  departure  of  the  locus  from 
the  pole.     Since      P         is  on  the  locus  the  sum  of  angles  must  apply. 

e4  +  e2  -  (e3  +  e    +  o5  +  e6  )    =    °  (4.69) 

solving  for 

96      =      f   e4    +   62    "    <e3    +   e4    +   66>    ]  (4«  70) 


As  point      P        approaches  the  pole,   the  angles  to  the  right  of  the 
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above  equation  become  simply  the  angle  subtended  between  the  various 

poles  and  zeros  and  the  complex  pole  in  question.     Obviously,   when 

the  angle  of  arrival  at  a  zero  is  determined  a  change  of  sign  is  required. 

4.5      Closed- Loop  Root  Locus 

Application  of  the  closed- loop  root-locus  assumes  that  the  poles 
and  zeros  of  the  open- loop  function  are  known  and  that  the  poles  of  the 
closed-loop  function  are  to  be  determined.     For  a  system  such  as  Fig- 
ure   4.  1    the  open-loop  and  closed- loop  function  are  related  by 

w<*>    ■     1    °G*1\  <4-71' 

Hence  the  poles  or  roots  of      W(z)      are  the  solution  of  the  equation 

1    +G(z)      =      0  (4.72) 

As    discussed  in  the  previous  Section,   both  functions  are  rational  poly- 
nomials in    z  .     Further,   the  zeros  and  gain  factors  of  both  functions 
are  identical.     Thus,   the  poles  of  the  closed-loop  function  are  the 
values  of    z    for  which 

/-  /    \  j  h  it 

G(z)      =      eJ 

h         =      1,     3,    5    ,  (4.73) 

The  precise  position  of  the  poles  or  the  root  locus  requires  that 

I  G(z)    |      =      1  (4.  74) 

As  in  the  case  of  the  open-loop  root  locus,   the  first  step  in  apply- 
ing this  graphical  root  finding  method  is  to  determine  the  locus  of  roots. 
The  approximate  root  locus  are  easily  found  by  application  of  the  rules 
which  have  been  presented  (6)  (15)    in  the  literature.     A  typical  closed- 
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loop  root  locus  is  shown  in  Figure    4.  15. 

4.  6      Additional  Design  Aids 

It  has  been  tacitly  assumed,   before  the  techniques  of  this    disser- 
tation are  employed,   the  optimum  system  performance  by  simple  gain 
adjustments  was  found  short  of  the  required  specifications. 

This  preliminary  gain  adjustment  can  be  aided  by  Figures    4.  16 

and   4.17.       These  Figures  relate  the  transient  response  quantities  of 

M      and      n        for  second  order  systems  with  no  position  error  to  the 
P 

location  of  the  closed-loop  complex  poles.     Figure    4.  16   applies  to  a 
system  with  sampler  only,   while  Figure    4.  17   applies  to  a  system  incor- 
porating a  zero  order  hold  network   (1)    (12)   .       Only  the  first  quadrant 
is  shown  since  complex  pole- zero  patterns  are  symmetrical  about  the 
real  axis  and  poles  in  the  second  quadrant  seldom  give  rise  to  a  satis- 
factory response.     The  figures  are  entered  on  the  circle  correspond- 
ing to  the  product  of  s -plane  open- loop  pole  and  sampling  period.     These 
circles  are  actually  the  closed- loop  root  locus  of  the  complex  poles 
and  proceeding  counter  clockwise  corresponds  to  increased  gain.     On 

the  circles  are  indicated  the  values  of   ,  M      and      n        that  a  second 

P 

order  system  with  a  pole  at  that  location  will  produce.  For  complex 
systems,  the  figures  can  be  used  on  the  basis  of  dominance  to  deter- 
mine if  the  various  values  of      M      and      n        will  satisfy  the  specifica- 

P 

tions  with    sufficient  margin  to  provide  for  the  deleterious  effect  of  the 
neglected  poles  and  zeros.     If  this  appears  possible,  a  more  exact 
study  can  be  made  using  the  closed-loop  root  locus.     The  figures  may 
also  be  helpful  for  making  the  first  selection  of  the  predominate  poles 
in  synthesis  process. 
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z  -plane 


pen -loop  P0^ 


Figure   4.  15.      A  typical  closed- loop  root  locus 
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Figure   4.  16.       Values  of    M    and    n      for  second  order  systems 


93 


S  plane,    pole  / 
Sgmphnq  Period 


Figure   4.17.     Values  of    M    and    n      for  second  order  system 
with  hold  P 
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The  dipole  required  to  raise  the  steady- state  error  constant  may 
require  such  close  spacing  as  to  render  praphical  methods  impractical. 
Applying  the  procedures  of  this  thesis  analytically  would  require  the 
finding  of  roots  of  polynomials.     Since  the  approximate  values  of  the 
roots  are  already  known,   this  is  not  as  difficult  a  task  as  may  first 
appear.     However,   an  alternate  method  avoiding  these  difficulties  works 
directly  with  the  open- loop  function.     A  pole  and  a  zero  are  introduced 
into  the  open- loop  function  of  proper  relative  values  to  give  the  desired 
increase  of  steady-state  error  constant  as  indicated  by  Equations    2.64, 
2.  65    and   2.  66.       The  expression  for  the  output  transform  can  be  found 
as 


C(z)      = 


z  k  A(z) 

"z~TT    X       k  A(z)   +   S(z) 


(4.  75) 


where      A(z)      and      B(z)      are  the  numerator  and  denominator  of  the 
open- loop  function  while    k     is  the  gain.     The  actual  output  at  the  sam- 
pling instants  are  readily  found  by  division. 

Unfortunately,   this  method  leaves  the  designer  in  doubt  as  to 
whether  the  final  output  time  response  will  be  satisfactory.     Actually 
the  above  procedure  is  not  so  involved  that  trial  and  error  is  out  of 
order.     Furthe  r,   as  an  aid  in  the  initial  selection  it  can  be  shown  that 
the  closer  the  loop  pole  of  the  dipole  is  to  unity  the  smaller  will  be  the 
residue  of  the  output  transfer  function  in  the  pole. 

This  effect  of  open- loop  pole  location  on  transient  response  when 
employed  to  increase  the  steady-state  velocity  error  coefficient  will 
be  illustrated  as  follows.     For  simplification  let  the    original  value  of 
K*        and  the  sampling  period  be  both  equal  to  unity.     Let  it  be  required 
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that  the  value  of  K*  be  raised  by  a  factor  of  five.  If  the  open-loop 
pole  and  zero  are  designated  by  a  and  q  and  the  closed- loop  pole  by 
p         then  it  follows  directly  from  Equations    2.65    and   4.44   that 

5      =      -L-^-9  (4.76) 

and 

1 
5      =       . . (4.77) 

1  +~ ,»-i — 

1-p  1-q 

Solving  Equation   4.  77   for    p    yields 

p  ■     4ttt-  <4-78> 

The  introduction  of  a  dipole  into  the  closed-loop  function  affects 
the  response  in  two  ways. 

(1)  The  original  terms  are  modified,  but  slightly  since,  as  a 
result  of  the  fraction  (z-q)/(z-p)  which  is  essentially  unity  at  the 
original  poles,   the  residues  at  these  poles  are  almost  unchanged. 

(2)  The  term  ,   _  £  ^ W  .       pn  is  added  . 

Consideration  of  this  term  will  show  that  with  movements  of  the 
dipole  near  the  point    1.0   the  coefficient  of  this  term  will  change  by 
a  factor  almost  directly  equal  to      (q-p)/(p-l)      .      Figure   4.  18, 
employing  the  previous  expressions,    gives  the  value  of  this  ratio  as 
the  open-loop  pole  is  moved  closer  to  unity  in  value. 
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Open- loop  pole    a 

Figure  4.  18.       The  effect  of  the  position  of  the  open- loop  pole 
on  the  output 
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CHAPTER   V 
APPLICATION  OF  SYNTHESIS  PROCEDURE 

5.  1       Introduction 

Chapters    III  and   IV   are  concerned  chiefly  with  presenting  the 
synthesis  philosophy  and  the  techniques  necessary  to  implement  the 
philosophy.     In  this  chapter  the  actual  process  of  synthesis  is  illus- 
trated by  a  variety  of  different  specifications.     Also  discussed  in 
detail  is  how  the  methods  of  this  dissertation  may  be  applied  to  sys- 
tems which  do  not  have  a  sampler  in  the  error  channel.     The  chapter 
is  concluded  by  a  discussion  on  synthesis  of  systems  with  multiple 
inputs. 

5.2      Complete  Example  of  Synthesis  Procedure 

This  section  presents  an  example  of  synthesis  using  the  method 
of  this  dissertation,  an  example  purposely  chosen  to  illustrate  the 
basic  elements  of  the  procedure.     Synthesis  is  initiated  by  determin- 
ation of  the  specifications  which  include  the  fixed  part  of  the  system 
and  the  time  domain  performance  characteristics.     The  specifications 
are  as  follows: 

(1)      The  plant  has  a  sampler  in  the  error  channel  operating  with 
period  of   0.  1    second  and  a  transfer  function  given  by: 


G    (s)      =  K    <s  +  19> (5.1) 

p  s(s+3.  7)(s+8)(s+25) 

(2)      The  time  domain  specifications: 
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Approximating  the  plant  as  a  second  order  system  and  employ- 
ing Figure  4.  1  5   as  explained  in  Section   4.  6    shows  that  specifications 
cannot  be  satisfied  by  adjustment  of  gain  and  a  compensator  will  be 
necessary.     Consequently,   the  synthesis  procedure  outlined  in  Figure 
3.  3  will  be  required. 

Two  items  which  must  be  decided  before   selection  of  the  closed- 
loop  pole-zero  configuration  are  the  difference  between  the  number  of 
z-plane  poles  and  zeros  of  the  plant  and  how  many  of  the  poles  and  zeros 
of  the  plant  are  to  be  retained  in  the  desired  transfer  function.     For 
the  first  item  it  is  seen  that  the  plant  has  one  sampler,    no  delays  and 
a  numerator  three  degrees  less  than  the  denominator,   hence  the    z- 
transfer  function  of  the  plant  will  have  one  more  pole  than  zero.     For 
the  second  item,    the  retention  of  plant  poles  and  zeros  is  largely  tied 
up  with  consideration  of  the  physical  equipment.     In  this  case  an  infin- 
ite value  for      K*        requires  the  pole  at  the  origin  to  be  retained  and 
arbitrarily  the  pole  at    3.  7   will  also  be  retained. 

With  these  constraints  the  simplest  pole-zero  configuration  can 

be  a  pair  of  complex  poles  and  one  real  zero.     Usirfg  Figure   4.  15 

again  as  a  guide  the  first  choice  of  pole- zero  selection  is  made  with 

poles  at      0.5     +  j  0.45        and  a  zero  at      0.3.     Using  the  correlation 

theorems 

M  .  21 

n  =        .  3 

P 
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n  =  8 

s 

K*        =         13 

v 

The  transient  response  requirements  are  met  with  sufficient  mar- 
gin to  permit  realization  of  a  higher  velocity  error  constant  by  intro- 
duction of  a  dipole.     First,   however,    the  open-loop  poles  must  be  ascer- 
tained.    This  is  readily  done  by  the  open-loop  root  locus  shown  in  Fig- 
ure   5.  1.       The  gain  is  fixed  by  including  one  pole  at  point    1.  0  and  then 
the  remaining  pole  is  found  at    0.  64.       Also  shown  in  this  figure  are  the 

construction  lines  used  for  measurements  and  angles  to  determine      n 

P 

and      M    .         n        is  estimated  by  reference  to  Figure    4.  7    . 
s 

The  s -plane  pole  of  3.  7   which  is  to  be  included  in  the  final  open- 
loop  function  corresponds  to  a  z-plane  pole  of    0.  691.       Since  the  first 
try  has  produced  a  pole  not  too  far  removed  from  this  location,   a 
closed- loop  root  locus  is  now  formed  based  on  the  zero  at    0.  3,    and 
poles  at    1    and     0.691    .       This  locus  is  shown  in  Figure    5.2    .     Choos- 
ing on  the  locus  poles  at        0.55   +  j  0.46        gives  the  following  time 
specifications 


M 

= 

.28 

n 
P 

= 

3 

n 

s 

s 

9 

K* 

V 

■ 

13 

With  the  gain  determined,   the  open-loop  function  is 


G(z)      =  .591  (z  -.3)  (5.2) 


(z  -    l)(z  -  .  691) 


In  accordance  with  the  discussion  in  Section  4.  6   an  open- loop 

pole  is  added  at      0.99      and  a  zero  at      0.95.       The  value  of      K*        is 

v 
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z-plane 

x    closed-loop  pole 
D  open- loop  pole 
0    zero 


Figure    5.  1. 


First  selection  for  pole- zero  configuration  and 
open-loop  root  locus. 
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Figure    5.  2.       Second  selection  for  pole-zero  configuration  and 
closed-loop  locus. 
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then  given  by 


K*        =  .591  (l-.95)(l-.3)      a      65>  (5>3) 

v  (1  -  .99)(1  -  .691) 


Further,    realization  of  the  open-loop  function  as  a  continuous  function 
without  time  delays   requires  a  zero  at  the  origin.     To  retain  the  dif- 
ference of  one  between  the  number  of  poles  and  zeros  an  additional 
pole  is  added  at      0.  001.       This  pole  was  purposely  placed  very  near 
the  origin  to  cancel  the  effect  that  the  zero  may  have  on  the  transient 
response.      The  desired  open-loop  system  is  thus  described  by  the 
transfer  function 


G(z)    =    K    J±l*)      =  .591  z(z  -  .95)(z  -  .3)  {54) 

B(z)  (z-l)(z-.99)(z-.691)(z-.001) 


The  inclusion  of  the  additional  poles  and  zeros  will  result  in  a 
shift  of  the  complex  poles  and  a  change  in  the  transient  response.      The 
new  pole  positions  could  be  determined  and  the  correlation  theorems 
employed  once  again.     But  a  simpler  method  would  be  to  form  the  out- 
put   z-transform  function      C(z)      to  a  unit  step  directly  from      G(z)    as 

C(z)      =      -5_     x    K-^Zl (5.5) 

z  -    1  K  A(z)    -    B(z) 

By  simple  division  the  output  at  the  sampling  instant  is  in  this  case 

C(t)      =      .  591  u   (t-.  1)    +    1.087  u    (t-.2)  +  1.  33  u    (t-.  3)  (5.6) 

+  1.  34  u   (t-.4)  +  1.212  u    (t-.  5)  +  1.065  u    (t-.6) 
+  0.966  uQ(t-.  7)  +  0.  935  uQ(t-.  8)  +  .  955  u    (t-.9) 


+  .995  u   (t-1.0) 
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The  transient  specifications  are  then  seen  to  be 


M 
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P 
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9 

Converting  the    z-transfer  function  to  its  equivalent  Laplace  trans- 
form with  poles  in  the  primary  strip  is  the  next  step.     This  involves 
taking  a  partial-fraction  expansion  of 

G(z) 
z 

and  multiplying  through  by     z     ,     converting  term  by  term  to  Laplace 

transform  with  the  aid  of  Table    2.  2      ,   and,    lastly,    collecting  the  terms 

over  a  common  denominator.      These  three  steps  take  the  form  shown 

below. 

~,    x            6.698z             5.514z              .9382z             .  246z  /c    ,. 

G(z)    =    -    -     -    (5.7) 

z  -    1  z  -   .99  z  -   .691  z  -  . 001 


nl    .          6.698             5.514            .938                 .246  ,,    _. 

G(s)     -   - (3.o) 

s  s  +  .l  s  +  3.  7  s-6.  9 


G(s)      =  20.95(s  +  .512)(s  +  15.94)  {^9) 

s(s  +  .  l)(s  +  3.  7)(s  +  69) 


The  compensation  is  now  determined  by  dividing  the  derived 
transfer  function,   Equation    5.9,    by  the  plant  Equation    5.  1,   yielding 

G'    (s)     =      2Q.95(s  +  .512)(s+8)(s  +  15.  94)(s  +  25)  (5 

(s  +  .  l)(s  +  19)(s  +  69) 
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The  network  is  not  realizable  in  its  present  form  but  the  addition 
of  a  fourth  pole  far  out  on  the  negative  axis,   has  a  negligible  effect  on 
the  transient  response,   and  yields  a  realizable  compensation  function 

G   (s)    =      20-  95(s+.  512)(s  +  8)(s  +  15.94)(s+25)  (5    n) 

c  (s  +  .  l)(s  +  19)(s  +  69)(s  +  100) 

This  function  is  realizable  in  terms  of  resistances  and  capacitors  by- 
standard  network  synthesis  procedures    (12)  (23)  . 

The  preceding  paragraphs  have  attempted  to  illustrate  the  basic 
synthesis  procedures.     However,    experience,    personal  preferences 
i    and  considerations  of  particular  systems  may  indicate  deviations  from 
the  initial  steps  outlined  above.     As  for  example,    if  the  entire  plant 
must  be  retained  as  part  of  the  final  desired  open-loop  function  it  may 
be  expedient  to  begin  the  design  on  that  basis.     This  would  correspond 
to  entering  the  Flow  chart  (Figure    3.  3)  initially  from  the  plant  posi- 
tion.     This  versatility  of  design  procedure  is  the  strong  point  of  this 
synthesis  method.     This  versatility  is  furnished  by  the  two  root  locus 
techniques  which  relate  the  movements  of  poles  in  one  function  on  the 
other  and  the  correlation  theorems  which  permit  interpretations  of  the 
pole-zero  configuration  in  terms  of  time   response. 

5.3      Further  Examples  of  System  Design 

In  this  section  three  systems  are  considered,    each  with  the  same 
continuous  fixed  components  involving  one  integration  and  an    s-plane 
pole  at     0.431.       The  first  system  includes  only  the  required  sampler, 
the  second  employs  a  sampler  and  holding  network,    while  the  third  is 
characterized  by  a  pure  time  delay  of  one  half  second,    giving  rise  to 
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the  following  three  transfer  functions, 


G,      (s)      = (5.12) 


IP 


s(s  +  .  431) 


1   -  e"S  K 


G7     (s)      =       — - x     " (5.13) 


"2p 


s(s  +  .431) 


G,     (s)      = x      e"   *5  S  (5.14) 


3p 


s(s  +  .  431) 


In  each  case  the   sampling  period  is  one  second. 

Besides   requiring  that  the  above  fixed  elements  be  included  in 
the  final  open-loop  transfer  function  each  system  must  also  satisfy 
the  time  specifications. 

M  <.        0.  5 

n  <         3 

P 

K  >        1 

v 

With  all  systems   required  to  meet  the  same  time  specifications 
at  the  sampling  instants  and  having  the  same  principle  s -plane  pole, 
the  same  closed-loop  z-transform  should  be  satisfactory  for  all  three 
systems.     A  pole-zero  configuration  with  one  pair  of  poles  at         0.  5  _ 
jO.  55      and  a  zero  at    0.  15    meets  the  above  specifications  w^th 
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Forming  the  output  transform  to  a  step  function  and  dividing  out 
gives  the  output  time  sequence  as 
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Forming  the  output  transform  to  a  step  function  and  dividing  out 
gives  the  output  time  sequence  as 
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c(t)      =      0.  65  u    (t-1)  +   1.  202  u    (t-2)  +  1.408  u     (t-3) 

+  1.306  u    (t-4)  +  .    .    .  (5.  15) 

* 
The  reason  for  the  difference  between  the  predicted  value  of      M      and 

its  actual  value  is  due  to  the  noninteger  value  for    n    as  explained  in 
Section   4.  2  . 

The  closed-loop  transfer  function  has  a  corresponding  open- 
loop  transfer  function 

=      .65(z  -   .15) (5>16) 

d  (z  -   l)(z  -   .65) 

It  should  be  noted  that  the  plant  pole  has  been  accommodated 

since 

-  .431  ,, 

e  =         .  65 

Although  to  this  point  the  procedure  has  been  the  same  for  all 
three  cases,   the  conversion  to  the   Laplace  transfer  function,    the  next 
step  in  deriving  the  compensating  network,    is  different  for  each  case 
and  will  be  illustrated  in  turn. 

(1)  For  the  first  system  the  plant  is  a  continuous  network  whose 
z-transform  is  characterized  by  a  zero  at  the  origin.  Adding  this  zero 
to  Equation  5.  15  and  a  pole  0.  001  to  maintain  the  proper  difference 
between  poles  and  zeros,    the  resulting  transfer  function  becomes 


Gld(z)      *      .65z(z-.15) (5.17) 

(z  -    l)(z  -   .65)(z  -   .001) 


Employing  the  methods  described  in  the  previous   section  yields 
the  equivalent  Laplace  transfer  function 
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Glrq(s)      =      1.6(s  +  2.82)  (5.18) 

s(s  +  .431)(s  +  6.9) 


the  transfer  function  of  the  compensator  is  then 


m  0.408(S  +  2.82)  {Sl9) 

d  (s    +   6.9) 


(2)      For  a  system  including  a  zero  order  hold  network  the  trans 
fer  function  of  the  open-loop  system  is  the  general  form. 


z    "    1       K      G      G      (z)  (5.20) 

so 


The         (z  -    1)/  z         term  corresponds  to  the  transfer  function  of 

-   sT 
1    -    e  ,     while  the  transfer  function      G    (s)      of  the  pure  integra- 

s 

tion  supplied  by  the  hold  network  and  the  transfer  function      G    (s)      of 
the  continuous  portion  of  the   system  result  in  the  z-transform 
G     G*?f)  (s)      in  the  form  of  Equation    5.  20    yields 


G,  .  (s)      =        Z  "    1       x zU  ?"    •  15) (5.21) 

a  z  (z   -   1)     (z  -  .65) 


and  then 

Gg  G20  (z)      =         Z(Z    "    *  15) (5.22) 

(z  -    l)2  (z  -   .65) 

The  conversion  of        G     G->^  (z )         to  its  corresponding   Laplace  trans 
form  is  similar  to  the  previous  example  yielding 
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hence 


r    i    \  r        i    \                 .449(s   +  1.514)  ,,    ?-. 

G    (s)  G?n  (s)      =      — ? S (5.23) 

20  sZ(s  +  .431) 


G,0(.)      =  8'98(S    +    1'514)  (5.24) 

6V  s(s  +  .431)(s  +  20) 


Where  the  pole  at    20.  0    has  been  added  to  enable  physical  reali- 
zation of  the  compensator 


(s)      =        0.0757(s    +    1.514)  (5>25) 

^C  (s    +    20) 


(3)      The  third  system  includes  a  pure  time  delay  of  one-half  the 
sampling  period.     Because  of  this  delay  a  zero  at  the  origin  is  not 
required  in  the  transfer  function.     Proceeding  directly  to  the  partial 
fraction  expansion  of  Equation    5.  16    gives 


Gd(z)      =    -i^ +      '— (5.26) 

z  -    1  z  -   . 65 


The  conversion  of  each  term  on  the  right  hand  side  of  the  above 
equation  to  the  s-domain  must  be  done  through  use  of  Table    2.  3      , 
Modified  z-transforms  ,    with    m    equal  to  one-half,   giving  rise  to 


G3d<s>  1.58  1.15  (5_2?) 


-  .  5  s  s  s    +   .431 

e 


Combining  terms  under  a  common  denominator,  adding  a  pole  at  twenty 
and  making  the  necessary  gain  adjustment  gives 
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m  8.58(s  +  1.52)  e"-5s  (5>28) 

*a  s(s  +  .431)(s  +  20) 


the  required  compensator  is 


(s)      .  0.075(s+1.52)  {5Z9) 

5C  (s    +   20) 


A  second  series  of  three  applications  evolves  when  the  value  for 
K*        is  increased  to  six  while  holding  all  other  specifications  the  same, 
The  open- loop  z-transfer  function  satisfying  these  conditions  is 


G„(«)        =     .65(z-.92)(2-.15)  (5>30) 

d  (z  -   l)(z  -  .  98)(z  -  .65) 


Employing  previously  discussed  methods  of  this  chapter  the 
required  open-loop  systems  are  described  by  the  transfer  functions. 


G        (g)  _        1.  574(s  +  .  083)(s  +  2.88)  (5.31) 

ld  s(s  +  .02)(s  +  .431)(s  +  6.9) 


G       *.)      =       l~e-S    x_9.477(s  +  .8)(s  +  1.43) ^  ^ 

a  s  s(s  +  .02)(s  +  .431)(s  +  20) 


G-,<.)         =        8.45(s  +  .083)(s  +  1.572)      /  -  5  s  (g>  33) 

s(s  +  .02)(s  +  .431)(s  +  20) 


The  output  time  sequence  to  a  unit  step  input  for  the  above  sys- 
tem is  found  to  be 
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c(t)    =    0.  65u    (t-1)    +    1.24uQ(t-2)     +     1.48uQ(t-3) 

+      1.  37  u   (t-4)    +    .     .     .  (5.34) 

o 


Each  of  the  six  systems  synthesized  in  this   section  were  instru- 
mented on  an  analog  computer.     The  computer  outputs  are  given  in 
Figures    5.  3   through    5.  8.     It  is  to  be  noted  that  the  output  at  sampling 
instant  checks  very  closely  with  the  calculated  output.     TJSble    5.  1  com- 
pares the  calculated  and  measured  values  for      K* 
r  v 

Liberty  was  taken  in  Figures    5.  5   and   5.  8   where  the  curves  are 
actually  the  input  to  the  delay  network  displaced  by  half  a  sampling 
period.      The  imperfection  of  the  computer  time  delay  required  this 
modification  of  the  actual  data.     Justification  for  this  procedure  as 
well  as  detailed  description  of  the  computer  study  is  given  in  Chapter 
VI. 

5.4      Sampled  Data  System  without  a  Sampler  in  the  Error  Channel 
The  sampled  data  systems  to  be  discussed  in  this   section  are 
characterized  by  a  lack  of  sampling  in  the  error  channel.     Figure    5.9 
is  the  block  diagram  of  such  a  system.     Reference  to  Table    2.2    gives 
as  an  expression  for  the  output 

R  G.    (z)  Q     (z) 

C(z)      =      i £ (5.35) 

1  +  H  Gl  Gz  (z) 

It  is  of  particular  interest  that  the  transform  of  the  input  signal  does 
not  appear  as  a  multiplicative  factor  in  the  numerator  of  the  expres- 
sion for  C(z)  .  Consequently,  it  is  not  possible,  as  was  done  pre- 
viously,  to  divide  out  the  transform  of  the  input  signal  and  form  a 

closed-loop  transfer  function  independent  of  the  input. 
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Table    5.  1.       Comparison  of  calculated  and  measured 
velocity  error  constants 


ystem 

Calculate 

id 

M 

easured 

sampler 

1.58 

1.68 

hold 

1.58 

1.64 

delay 

1.58 

1.  55 

sampler 

6.2 

5.4 

hold 

6.2 

5.7 

delay 

6.2 

6.0 
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Figure   5.9.       A  sample-data  control  system 
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By  the  expediency  of  dividing  both  sides  of  Equation    5.  35   by  a 
z-transform  of  the  input  a  closed-loop  transfer  function  can  be  defined 
as 

r  R  G.    (z)  Q     (z)   /  R(z) 

(z)        =      i L- (5.36) 


R 


1    +    H  Gl  G2  (z) 


It  should  be  noted  that  the  transfer  function  applies  only  for  a  specific 
input  and  for  this  reason  will  be  referred  to  as  a  limited  closed-loop 
transfer  function  with  designation       W     (z)    . 

In  general,   limited  closed-loop  transfer  functions  do  not  have  the 
same  zeros  as  the  open-loop  transfer  functions  even  though  the  system 
provides  a  unity  feedback  path.     Because  of  this  characteristic  it  is 
not  possible  to  employ  the  open- loop  root  locus  for  finding  the  poles  of 
the  open-loop  functions  nor  are  the  expressions  for  the  steady  state 
error  constants  in  terms  of  closed-loop  poles  and  zeros  applicable. 
However,  the  response  to  any  particular  input  can  be  evaluated  dir- 
ectly in  the  manner  employed  in  Section   4.  2   for  a  unit  step  and  the 
response  shaped  as  desired  by  proper  selection  of  the  closed-loop  poles 
and  zeros.     The  steady  state  error  can  be  found  by  direct  application 
of  the  final  value  theorem,  Equation   2.  55.     For  systems  with  zero 
steady  state  error  to  a  step  input,   Equations    4.  14   and   4.  19   are 
directly  applicable  for  synthesis  to  a  prescribed  transient  response. 

With  the  limited  closed-loop  transfer  function  determined,  the 
desired  open-loop  transfer  function  may  be  found  by  a  method  such  as 
the  following.       With  a  step  input  and  setting      H(s)      to  unity  the  remain- 
ing expressions  defined  in  Equation   5.  36  are  polynomials  in      z      with 
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expressions 


R(z) 


(5.37) 


RG1  (z)      = 

krl  arl   (z) 

(z-1)  b1   (z) 

fl       (y)           = 

k2a2  (z) 

u2  [    ' 

b2(z) 

Gl  G2  (z)    = 

k12a12  {z) 

bj   (z)  b1  (z) 

(5.38) 


(5.39) 


(5.40) 


R 


U) 


K  Q  (z) 
P  (z) 


(5.41) 


In  Equations    5.38   and    5.40    ,   which  are  z-transforms  of  the  product 
of  two  Laplace  transforms,   the  convenient  property  that  s-plane  poles 
have    a   one   to   one    image  on  the    z-plane  is  utilized  by  factoring  the 
denominator  into  the  parts  corresponding  to  each  Laplace  transform. 
Substitution  of  the  above  expression  into  Equation   5.36    yields 


K 


QU) 


krl  arl(z)  k2a2(z) 


b^z)  b2(z)     +    k12  a12  (z) 


(5.42) 


Generally,    either      G.    (s)      or      G?  (s)      are  known  in  part  or 
whole.     Assuming  that      G.    (s)      is  fixed,   then  the  zeros  of      G?  (s) 
are  fixed  by  Equation   5.42    ,  after  the  closed- loop  zrros  and  gain  con- 
stant are  selected.       Consequently,   the  factors  of     b?   (z)      must  be 
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chosen  to  satisfy  Fquation   5.  39   and  also  combine  with  the  remaining 
open-loop  poles  and  the  open-loop  zeros  to  form  the  poles  of  the  closed- 
loop  function  with  the  exception  of  the  pole  at  the  origin  introduced  by 
the  transform  of  the  input. 

A  first  selection  of  the  poles  is  aided  by  the  consideration  of  the 
physical  system  and  the  type  of  compensation  to  be  employed.     A  phase 
lead  network  requires  a  pole  near  the  origin.     A  dipole  for  reduction 
of  the  steady  state  error  requires  a  pole  near  the  point      1.0    .     The 
total  number  of  open-loop  poles  is  one  less  than  the  number  of  closed- 
loop  poles.       The  expression  for      A?(z)      can  be  solved  for  by  equating 
the  numerators  of  each  side  of  Equation    5.  42   .     A  first  trial  expres- 
sion for      G2  (z)      is  then  known  and  function      G~  (s)      can  be  deter- 
mined by  methods  illustrated  in  the  previous  sections  of  this  chapter. 
Then  the  expression  for       G.  G?  (z)       may  be  formed  and  by  root  locus 
technique  the  closed-loop  poles  found.       In  general,   these  poles  will 
not  be  those  previously  selected.       However,   the  resulting  poles  may 
very  well  satisfy  the  specifications  thereby  completing  the  design. 
If  this  is  not  the  case  the  process  must  be  iterated.       The  precise 
method  for  determining        G,  (s)        may  vary  with  different  system 
arrangements,   but  with  experience  systems  of  this  sort  can  be  readily 
handled  by  application  of  the  above  principles. 

5.  5      Systems  with  Multiple  Inputs 

Because  of  the  particular  application  or  because  of  disturbances, 
all  feedback  control  systems  have  multiple  inputs.     For  example.    Fig- 
ure   5.  10    shows  the  block  diagram  of  a  sampled  data  servo  system  in 
which  the  output  is  corrupted  by  a  spurious  signal      u(t)      entering  the 
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Figure    5.  10.       Tandem  single- loop  system  with  secondary 
input 
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system  at  point    A  .       Thus  a  feedback  control  system  is  characterized 
by    n    closed- loop  transfer  functions,    the  customary  or  primary  trans- 
fer function  relating  output  with  the  reference  input  and  a  number  of 
secondary  transfer  functions  for  each  secondary  input  to  reference  out- 
put.      Whatever  the  source  of  the  secondary  input  or  wherever  the 
point  of  entry  to  the  system,    satisfactory  design  involves  the  realiza- 
tion of  suitable  dynamic  characteristics  for  each  of  the  possible  trans- 
missions. 

In  Figure    5.10       G?  (s)      and      G.,   (s)      may  represent  transfer 
functions  of  fixed  parts  of  the  system  while      G     (s)      is  the  transfer 
function  of  the  compensator.       For  a  signal      u(t)      with  transform      U(z) 
the  expression  for  the  transform  of  the  output  is 


UG     (z) 
C(z)      =    ^ (5.43) 

1    +  Gc  G2  G3  (z) 


All  factors  on  the  right  side  of  the  above  Equation  have  been 
fixed  by  synthesis  of  the  primary  transfer  function  given  in  the  speci- 
fications.      No  provision  is  available  for  modifying      C(z)      if  the  results 
of  Equation   5.  43   are  unsatisfactory.       More  latitude  is  offered  if  the 
compensator  transfer  function      G     (z)      is  split  into  tandem  and  shunt 
sections.       This  arrangement  is  illustrated  by  the  system  in  Figure 
5.  11a   .       Figure    5.  lib    shows  the  same  system  reduced   to  a  single 
loop.     The  transfer  functions  of  the  blocks  are  related  by  the  equations 


Ga  (s)      =       Gj   (s)  G2  (s)  (5.44) 
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(a)      System  with  two  input  employing  feedback  and  tandem  compensator 


R(s)  y^  E(s)^ 


Ga  (*> 


Gb  (s) 


(b)      Tandem  single- loop  system 


Figure    5.11.       Two  equivalent  systems 


C(s) 
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G,  (s) 

G     (s)      =      ± (5.45) 

1  +  G2  (s)  G     (s)  G     (s) 


The  closed-loop  z-transform  transfer  function  for  the  system  in  Fig- 
ure   5.11b    is  given  >by 


C                               GaGb(z) 
■i-  (z)      =     a— 5 u  (5.46) 

*  1    +  G     G,    (z) 

a      b  v 


Figure    5.  12a    shows  the  system  of  the  previous  figure  redrawn 
and  the  reference  input  omitted  to  emphasize  the  transfer  character- 
istics of  the  system  between  secondary  input  and  reference  output.     By 
block  diagram  manipulation  the  system  can  be  reduced  to  that  shown 
in  Figure    5.12b     where  the  transfer  functions      G     (s)      and       G,    (s) 

3.  D 

are  the  same  as  defined  by  Equations    5.44   and   5.45   .     The  transfer 
function  of  the  output  is  given  by  the  expression 

U  G     (z) 

C(z)      =       2 (5.47) 

1    +   G     G,    (z) 
a     b 

Using  the  methods  of  the  previous  section  the  limited  closed- 
loop  transfer  function  for  the  secondary  input  is  obtained  as  follows 


c                        U  G.    (z)   /  U(z) 
■ijr   (z)       =    2 (5.48) 

1    +   G     G,    (z) 

a      b  N    ' 


Noting  that  the  denominators  of  Equations    5.46   and    5.48  have  a  com- 
mon term  permits  as  a  more  convenient  expression  for 
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Figure    5.  12.       Two  equivalent  systems 
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u 


(z) 


U  Gb  (z) 


U(z)  Ga  Gb  (z) 


R 


(z) 


(5.49) 


Equation    5.49   places  in  evidence  the  possibility  of  satisfying 
independently  specifications  for  both  the  primary  and  secondary  trans- 
mission characteristics.     However,       U(z)      and      R(z)      must  be  the 

q 

same  function  if  the  equation  is  to  be  meaningful.     Poles  of  — =^ —  (z) 

Q 

and         — = —   (z)        are  identical  with  the  exception  of  an  additional  pole 

Q 

at  the  origin  of  the  expression  for        — « —  (z)        .     Consequently,   the 

Q 

design  freedom  for  meeting  specification  on         -f    (z)         i-s  limited  to 
choice  of  gain  constant  and  zeros  of        U  G,    (z)    . 

The  choice  of  zeros  and  gain  constant  for  meeting  specifications 
are  aided  by  the  relationships  given  in  Section   4.  2    .     With  the  design 


of  the  system  for 
C 


R 


(z )        already  completed  and  with  the  zeros 


of 


U 


(z)         selected  Equation    5.  49    may  be  solved  for  the  remain- 


ing unknown        U  G,    (z)      .     Hence 


U  Gb  (z)      = 


U 


(z) 


R 


(z) 


U(z)   Ga  Gb  (z) 


(5.50) 


By  procedures  described  in  detail  earlier  in  this  chapter  the  transfer 
function  for      G,    (s)      may  be  evaluated.     The  transfer  function  of  the 
feedback  compensator       G  .  (s)      is  then  given  by 


G4  (s)      = 


G3  (s) 
Gb(s) 


G3  (s)  G2  (s) 


(5.51) 
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It  is  apparent  that  the  poles  of       G.,  (s)       and       G^  (s)        ,   which  are 
also  poles  of  the  open- loop  function,  are  zeros  of       G.   (s)     .       Conse- 
quently      G.  (s)        will  in  general  require  one  or  more  poles  at  infinity. 
Hence,    in  practice       G.  (s)       can  only  be  realized  approximately  by      R- 
C    circuits  and  high  gain  amplifiers. 

The  expression  for  the  transfer  function       G.    (s)      ,   given  in 
terms  of  the  previously  determined  desired  open-loop  transfer  func- 
tion,       G,  (s)      ,     is 
a 

Gd(s) 

G.    (s)      =    2 (5.52) 

1  G2  (s)    G     (s) 


5.  6      Suggestions  for  Future   Work 

The  field  of  sampled-data  theory  and  application  is  still  in  its 
infancy.       Many  problems  remain  to  be  studied.       One  field  that  should 
prove  very  useful  is  the  use  of  sampled-data  theory  in  the  analysis  of 
continuous  linear  and  nonlinear  systems.     In  this  case  continuous  func- 
tions would  be  represented  by  sampled  approximations.     The  analysis 
of  sampled-data  systems  with  nonlinear  components  offers  ample  pos- 
sibilities for  study.       These  areas  of  sampled-data  theory  are  not  dis- 
cussed in  this  dissertation. 

Of  the  problems  which  are  a  direct  outgrowth  of  this  dissertation 
three  are  of  particular  interest:    first,    extention  of  this  synthesis  pro- 
cedure for  systems  with  non-periodic  operation  of  the  various  sam- 
plers in  the  systems  and  consideration  of  wide  sampling  pulses  which 
cannot  be  replaced  by  impulses.      Second,   further  extention  of  the 
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procedure  to  control  the  shape  of  the  response  between  the  sampling 
instants.     The  third  problem  concerns  the  need  to  study  the  zeros  of 
a  z-transform  to  develope  a  method  of  predicting  the  location  of  zeros 
from  the  time  function. 
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CHAPTER   VI 
DESCRIPTION  OF  THE  ANALOG  COMPUTER  SET-UP 

An  electronic  analog  computer  study  was  undertaken  to  verify  the 
synthesis  methods  of  this  dissertation.     The  computer  was  set  up  to 
correspond  on  a  transfer  function  basis  with  the  physical  system. 
This  is  shown  by  Comparing  Figure    6.  1,    where  a  block  diagram  of  a 
sample-data  position  servo  is  shown,   and  Figure    6.  2    ,  where  its  com- 
puter analog  is  shown.     The  system  consists  of  a  sampler,   phase  lead 
network  and  a  motor.     In  practice  the  latter  device  may  be  a  commu- 
tator or  any  of  a  variety  of  analog  to  digital  converters    (23)  (24).     A 
photograph  of  the  entire  computer  layout  is  shown  in  Figure    6.  3    .     The 
important  components  of  the  set-up  are  identified  on  the  figure. 

The  computer  set-up  on  a  transfer  function  representation  of  a 
sampled  data  system  served  as  more  than  just  an  aid  to  calculation, 
but  as  a  direct  analogy  of  the  physical  system  and  consequently  a  check 
on  the  analytical  methods  of  this  dissertation.       Since  the  methods  of 
this  thesis  are  based  on  linear  theory,   the  electronic  analog  computer 
simulation  provides  a  more  rigorous  test  than  a  physical  system.     No 
excuse  is  possible  that  nonlinearities  obscured  the  results  predicted 
by  the  theory.     Also,   the  versatility  of  the  computer  permitted  the 
testing  of  more  possible  system  variations. 

As  well  as  the  system  described  above,   a  second  system  using  a 
hold  network  and  a  third  system  incorporating  a  pure  time  delay  as 
well  as  the   sampler  were  tested.       Each  of  the  systems  were  designed 
to  meet  the  same  specifications  and  consequently  the  z-plane  pole-zero 
configuration  was  the  same.     Two  compensating  networks  were  provided 
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Figure    6.  1.       Simple  sample- data  servo  system 
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for  each  system--the  first  included  only  a  phase  lead  network  to  adjust 
for  rise  time  and  peak  overshoot  while  the  second  added  a  dipole  to 
increase  the  velocity  error  constant.      The  details  of  the  procedure  for 
design  of  the  compensators  were  explained  in  Section    5.  3    .     Also, 
shown  in  Figures    5.  3   through   5.  8   are  copies  of  the  output  voltage  ver- 
sus time  curves   secured  from  the  computer  set-up.       These  curves 
along  with  the  tabulation  of  values  for      K*        in  Table    5.  1    represent 
the  observed  data  of  the  computer  study.       The  very  close  check 
between  the  analytical  results  and  the  observed  results  verified  the 
methods  of  this  dissertation  and  no  further  consideration  of  this  mat- 
ter will  be  made  in  this  chapter.       Rather,   this  chapter  will  limit  its 
discussion  to  the  new  components  which  were  developed  to  simulate 
the  various  system  characteristics  which  are  peculiar  to  a  sample- 
data  control  system.       The  components  referred  to  are  (1)   the  Sam- 
pling unit,     (2)   the  holding  circuit,    (3)   the  pure  time  delay,   and    (4) 
the  input  signal  synchronizer  circuit. 

(1)      A  front  view  of  the  Sampling  unit   (16)    is  shown  on  the  left  of 
the  Figure    6.  3    .     Four  sets  of  terminals  are  shown  since  actually  this 
particular  device  furnished  four  sampling  switches.       A  top  view  clearly 
exhibiting  the  various  parts  making  up  the  device  is  available  in  Figure 
6.4    .     From  this  figure  it  is  seen  that  a  D.  C.    motor*  drives  four  ganged 
differential  gears  through  a  step  down  gear  train.     The  output  shaft  of 
each  of  the  four  differential  gears  drives  a  split  cam.     Adjustment  of 
the  cam  permits  control  over  the  switch  closure  time.       Ford  automo- 
bile distributer  contacts  are  used  for  the  actual  switch.     Each  differ- 


*Electric  Indicator  Co.  ,    1/25  HP.    115  volts,   style  ASC  motor. 
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ential  gear  unit  is  furnished  with  a  knob  which  permits  displacing  the 
angle  alignment  between  its  input  and  output  shafts.       Since  the  input 
of  each  differential  is  driven  off  the  previous  differential  output  (except 
for  the  first)    this  permits  relative  phasing  between  the  closures  of  the 
four  sampling  switches. 

The  Electronic  Counter*    shown  in  Figure    6.  3    was  used  to  mea- 
sure the  sampling  period  and  sampling  duration.       The  sampling  switch 
contacts  were  readily  transferable  from  the  computer  circuit  to  the 
input  of  the  counter  in  parallel  with  a  battery  and  resistor  circuit  by 
means  of  the  knife  switch  shown  in  front  of  the  Sampling  unit.     This 
permitted  control  of  the  sampling  period  by  adjustment  of  the  potenti- 
ometer in  series  with  the  D.  C.    motor  field.       The  electronic  counter 
was  adjusted  to  give  time  interval  measurements  to  a  tenth  of  a  milli- 
second.     A  one  second  sampling  period  was  chosen  and  it  was  found 
that  the  motor  speed  was  sufficiently  steady  to  hold  the  sampling  period 
variation  to  less  than   _  5   milliseconds  during  a  computer  run. 

The  closure  time  of  the  sampling  switch  was  set  for    59  milli- 
seconds.      Employing  Equation   2.  2    indicates  that  the  actual  gain  of 
the  computer  simulated  system  needed  to  be     16.95    times  the  calcu- 
lated gain  in  order  that  the  pulse  correspond  in  strength  to  an  impulse. 
That  the  system  would  react  to  this  pulse  in  approximately  the  same 
manner  as  to  an  impulse  is  assured  since  the  principle  time  constant 
of  the  simulated  system  was    2.  3    seconds.       The  sampling  period  of 
one  second  and  the    59    milliseconds  pulse  were  chosen  as  a  compro- 
mise between  the  gain  and  impulse  approximation  requirements. 


♦Hewlett  Packard    522B   Electronic  Counter. 
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(2)  The  sampler  plus  the  holding  circuit  was  simulated  by  the 
circuit  depicted  in  Figure    6.5.       It  has  the  advantage  of  employing 
only  one  computer  amplifier  and  the  switch  was  furnished  by  the  sam- 
pling unit.       For  proper  operation  the  circuit  functions  as  follows: 
when  the  switch  closes  the  amplifier  output  voltage  very  quickly 
assumes  the  same  value  as  the  input  (with  reversed  sign);  when  the 
switch  opens  the  output  voltage  due  to  action  of  the  capacitor  will 
remain  constant  until  the  switch  closes  next  time.       It  should  be  noted 
that  the  closure  time  of  the  switch  should  be  ample  to  insure  that  the 
capacitor  will  charge  up  to  its  correct  value  of  potential.     With  a 
switch  closure  time  of    15   milliseconds  suitable  parameter  values  for 
R      and      C      were  found  to  be      0.  01    megohms  and    0.015      microfa- 
rods  respectively.     The  preceeding  values  give  a  time  constant  of   0.  15 
milliseconds.       Since  the  sampling  switch  is  closed  for  ten  time  con- 
stants it  can  be  assumed  that  the  difference  between  the  input  and  out- 
put signal  magnitudes  is  negligible. 

(3)  Pure  time  delays  or  transportation  delays  can  only  be  approx- 
imated by  linear  continuous  networks.       For  this  computer  study  the 
Pade  fourth  order  approximant    (12)   was  simulated  by  the  computer 
circuit  shown  in  Figure    6.6    (25)    .        The  value  of  time  delay   X    was 
set  for    0.  5    seconds  and  the  six  operational  amplifiers  were  secured 
from  the  separate  Donner  unit  shown  on  the  right  of  Figure    6.  3    .     In 
Figure    6.  7  the  output  of  this  circuit  to  an  input  ramp  signal  is  shown. 

It  is  significant  that  the  output  is  not  identically  zero  during  the  time 
delay  but  rather  performs  a  low  amplitude  oscillation.      The  oscilla- 
tion is  more  pronounced  if  the  input  signal  contains  a  large  proportion 
of  high  frequencies.       For  this  reason,   the  delay  circuit  was  placed 
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Figure   6.  5.       The  sampler  plus  the  holding  circuit  are 
simulated  by  this  circuit 
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Figure   6.  7.       Input  and  output  signals  of  the  simulated 
time  delay 
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at  the  output  of  the  open- loop  system  permitting  the  system  time  con- 
stant to  reduce  the  magnitude  of  the  complementary  frequency  compon- 
ents introduced  by  the  sampler.       Nevertheless  the  output  of  the  delay 
still  contained  an  oscillation  not  present  in  the  input  signal.       However, 
the  system  operates  open  circuit  during  this  period  and  the  signals  do 
not  interfere  with  the  proper  operation  of  the  system  since  the  delay 
assumes  the  proper  output  value  by  the  time  the  sampling  switch  closes 
again. 

(4)      Calculation  of  transient  responses  by  analytical  methods  is 
usually  performed  with  both  the  initial  conditions  and  the  input  function 
zero  at  time  zero.       For  a  sampled  data  system  this  requires  input 
functions,    such  as  ramps,   to  be  first  applied  at  a  sampling  instant. 
For  a  computer  study  this  would  require  that  the  start  of  a  computer 
run  be  synchronized  with  a  sampling  instant.       The  circuit  depicted  in 
Figure    6.  8   was  found  to  provide  this   synchronization.       This  circuit 
utilizes  a  second  sampling  switch  A  from  the  Sampling  unit  which  must 
be  synchronized  with  the  switch  used  for  sampling  in  the  computer. 
The  start  and  stop  switch  are  both  closed  when  a  computer  run  is  to 
be  initiated.       Then  upon  closure  of  the  sampler  contacts     A   ,    relay 
contact    B    closes  a  series  circuit  consisting  of  the  coil  and  battery 
and  thus  locking  the  relay  in  the  operate  position  while  closure  of  relay 
contacts    C     starts  on  the  computer  run.       The  Donner  computer  fur- 
nishes external  operate  contacts  which  were  thus  employed.       The  com- 
puter run  is  terminated  by  placing  the  stop  switch  in  the  open  position. 
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Figure   6.  8.       Input  signal  synchronizer  circuit 
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